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Matroid Polytopes: Algorithms, Theory and AppHcations. 

Abstract 

This dissertation presents new results on three different themes all related to matroid 
polytopes. First we investigate properties of Ehrhart polynomials of matroid polytopes, 
independence matroid polytopes, and polymatroids. We prove that for fixed rank their 
Ehrhart polynomials are computable in polynomial time. The proof relies on the geometry 
of these polytopes as well as a new refined analysis of the evaluation of Todd polynomials. 

Second, we discuss theoretical results regarding the algebraic combinatorics of matroid 
polytopes. We discuss two conjectures about the /i*-vector and coefficients of Ehrhart poly- 
nomials of matroid polytopes and provide theoretical and computational evidence for their 
validity. We also explore a variant of White's conjecture which states that every matroid 
polytope has a regular unimodular triangulation. We provide extensive computational ev- 
idence supporting this new conjecture and propose a combinatorial condition on simplices 
sufficient for unimodularity. Lastly we discuss properties of two dimensional faces of matroid 
polytopes. 

Finally, motivated by recent work on algorithmic theory for non-linear and multicrite- 
ria matroid optimization, we have developed algorithms and heuristics aimed at practical 
solutions of large instances of these difficult problems. Our methods primarily use the lo- 
cal adjacency structure inherent in matroid polytopes to pivot to feasible solutions which 
may or may not be optimal. We also present a modified breadth-first-search heuristic that 
uses adjacency to enumerate a subset of feasible solutions. We present other heuristics, 
and provide computational evidence supporting these new techniques. We implemented all 
of our algorithms in the software package MOCHA (Matroids Optimization Combinatorial 
Heuristics and Algorithms). 



-iii- 



Acknowledgments 

First and foremost to my advisor Jesus De Loera who started me on this journey as an 
undergraduate. He guided me around the pitfalls, encouraged me, mentored me, provided 
me with problems and kept challenging me at every turn. His patience was remarkable and 
he was my number one cheerleader. Most of all, he taught me that mathematics without 
passion is not math at all. 

To my wife Tami Joy whose support and love made this all possible. Her utmost 
confidence in me is amazing. It lifted my spirits and kept me going. She can always make 
me laugh, and is always there to lean on. 

My thanks to Matthias Koppe for sharing his mathematical and computational skills. 
His careful and thoughtful approach to mathematics set a great example and he always 
made himself available for any math or computer questions. 

To my family and friends for all their support, understanding and love. My thanks to 
the UC Davis mathematics department staff who were very helpful and instrumental in my 
success. 

Thanks to Francisco Santos for all our insightful conversations on triangulations of 
matroid polytopes and more. 

My appreciation to Eric Babson for serving on my qualification exam and dissertation 
committee. 

Support for this dissertation was provided by VIGRE (DMS-01-35345, DMS-0636297), 
University of CaHfornia at Davis (DMS-0608785) and a gift from IBM. 



-iv- 



CHAPTER 1 

What are Matroids and who are their polytopes? 

1.1. Matroids 

Matroids are combinatorial objects that naturally encapsulate the idea of independence. 
They are found in matrices, graphs, transversals, point configurations, hyperplane arrange- 
ments, greedy optimization, and pseudosphere arrangements to name a few. One of the 
reasons matroids have become fundamental objects in pure and applied combinatorics are 
their many equivalent axiomatizations. For an excellent guide and thorough treatment of 
matroids we recommend [0x192, Wel76, Sch03], which we refer to for most all of our def- 
initions and theory. We now offer a brief introduction to matroids and their associated 
polytopes. We begin with the most basic and prominent matroid definition (axiomatiza- 
tion). First we adopt the notation X + y := XiJ {y} and X — y := X \ {y}. A pair {S,I) is 
called a matroid if is a finite set and T is a nonempty collection of subsets of S (denoted 
2'^) satisfying: 

{i) if / G T and J C /, then J G J 

(1.1) 

{a) I, J £ I and |/| < | J|, then I + z £ I for some z £ J\I. 

For a matroid Ai = {S,I), I S is independent if / G X and dependent otherwise. B C S 
is called a base or basis if it is an inclusion maximal independent set. That is, S G X and 
there does not exist A £ I such that B C A (1 S. We denote the collection of bases of M 
as or sometimes B when the matroid is clear. An easy consequence of Equation 1.1 is 
that all bases have the same cardinality: If not, let Bi, B2 be two bases of {S,2) such that 
\Bi\ < 1^2 1- Then by (ii) in Equation 1.1 there exists z £ B2 such that Bi + z £ I. But 
this contradicts the inclusionwise maximality of Bi . 

Given any [/ C S the rank of U is given by the rank function of A^, tp: 2["1 ^ Z where 

ip{A) := max{ \X\ \ X O A, X £l}. 
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The common cardinality of the bases of Ai is called the rank of the matroid. Throughout 
we will refer to the rank of M as rank(A^) := rank(S'). 

Example 1 (Uniform matroid). Let [n] := {1,...,?7-}. The uniform matroid of rank 
r is := {A CI [n] \ \A\ < r}. It is easy to see that [/'"'" satisfies both (?) and (ii) of 
Equation 1.1. 

i72'4 = {{l,2},{l,3},{l,4},{2,3},{2,4},{3,4}} 

We state an additional useful axiomatization on the independence sets. 

Lemma 2 (Theorem 39.1 [Sch03]). Let S be a finite set and let I be a nonempty collection 
of subsets satisfying Equation l.l{i). Then Equation l.l{ii) is equivalent to: 

(1.2) i/I, J gT and |/\ J| = 1, |J\/| = 2, then I + z £l for some z € J\I. 

Proof. Obviously Equation l.l(n) implies Equation 1.2. Conversely, Equation 1.1(m) 
follows from Equation 1.2 by induction on |/\ J|, the case |/\ J| =0 being trivial. If 
|-^\ «^| ^ 1) choose i £ I\J. We apply the induction hypothesis twice: first to I — i and J 
to find j £ J\I with I — i+ j G T, and then to I — i + j and J to find j' G J \ (/ + j) with 
I — i + j + j' £1. Then by Equation 1.2 applied to / and I — i + j + j', we have I + j £ I 
or I + f el. □ 

We can now prove one of the most useful axiomatizations for the purposes of this 
dissertation. Considering the properties of Equation 1.1 it is natural to think that the 
bases of M are all that are needed to fully describe a matroid and indeed the following 
theorem gives an important characterization of matroids in terms of bases. 

Theorem 3 (Theorem 39.6 [Sch03]). Let S be a set and B be a nonempty collection of 
subsets of S. Then the following are equivalent: 

(i) B is the collection of bases of a matroid; 

(1.3) (ii) if B,B' eB and xe B'\B, then B' - x + y e B for some yeB\B'; 
(Hi) if B, B' £ B and x £ B' \B, then B' — y + x £ B for some y £ B\ B' . 

Proof, (z) =^ (ii): Let B be the collection of bases of a matroid {S,I). Then all sets 
in B have the same size. Now, let B,B' £ B and x £ B' \ B. Since B' — x £ I, there exists 
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ayeB\B' with B" := B' - x + y £ 1. Since \B"\ = \B'\, we know B" G B. 

{Hi) =^ (i): {Hi) directly implies that no set in B is contained in another. Let X be the 
collection of sets / with I B for some B €z B. We check Equation 1.2. Let I, J £ I with 
|I\ J| = 1 and |J\/| = 2. Let / \ J = {x}. 

Consider sets B,B' £ B with I B, J CI B' . If x G B' , we are done. So assume x ^ B' . 
Then by {Hi), B' — y + x £ B for some y £ B'\B. As | J \ /| = 2, there is a z G J \ / with 
z ^ y. Then I + zCB' — y + x, and so / + z G X. 

{a) =^ (iiz): By the foregoing we know that {Hi) implies {ii). Now axioms {ii) and 
{Hi) interchange if we replace B by the collection of complements of sets in B. Hence the 
implication {ii) ^ {Hi) also holds. □ 

Parts {H) and {Hi) of Equation 1.3 might seem identical, but {ii) says that we can pick 
an element to remove from B' and there exist some element of B to add to get a basis while 
{Hi) says that we can pick an element to add to B' and there exists some element of B to 
remove to get a basis. Having the flexibility to choose the element to add or remove is a 
useful property for many proofs, algorithms and heuristics on the bases of matroids. 

Example 4 (Graphical Matroid). Let G = {V, E) be graph with nodes V and edges E. 
The collection of spanning forests of G are the bases of the graphical matroid on ground set 
V, denoted Mq- We can verify the basis exchange axiom Equation 1.3{iii) by considering 
two spanning forests Fi,F2. If e G -F2 \ -^i, then Fi + e contains a cycle C in G. Also there 
exists / G C \ F2 and thus Fi — / + e is a spanning forest. 

Matroids play a prominent role in combinatorial optimization and are intrinsically tied 
to the greedy algorithm which we state now: If it; : 5* ^ M is a weight function we use the 
notation w{I) := X^jg/ w{i) where ICS. Letting T C 2'^ be closed under taking subsets 
and given a weighting w, the greedy algorithm starts by setting I := 0, and repeatedly 
choosing y £ S\I such that I + y £ I where w{y) is as large as possible. We terminate if 
no such y exists. 
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Figure 1.1. (a) A graph G and (b) its spanning forests, i.e. bases Bmg -^G- 

For arbitrary X C 2'^ the greedy algorithm does not always terminate at the maximum 
solution, although for matroids this will always be the case. In fact, the converse is true: If 
the greedy algorithm always works on Z then it is a matroid. 

Theorem 5 (Theorem 40.1 [Sch03]). LetZ he a nonempty collection of subsets of a set 
S, closed under taking subsets. Then the pair (5*, X) is a matroid if and only if for each 
weight function w : S ^ M+, the greedy algorithm leads to a set I £ I of maximum weight 
w{I). 

Proof. Necessity. Let {S,2) be a matroid and let w : S ^ M+ be any weight function 
on S. Call an independent set / good if it is contained in a maximum-weight basis. It 
suffices to show that if / is good, and y is an element in S\I with I + y £ I and with w{y) 
as large as possible, then I + y is good. 

As / is good, there exists a maximum-weight basis B ^ I. y £ B, then / -|- y is 
good again. If y ^ B, then there exists a base B' containing I + y and contained in B + y. 
So B' = B — z + y for some z £ B\I. As w(y) is chosen maximum and as / -|- z £ Z 
since I + z Q B, we know w(y) > w{z). Hence w{B') > w{B), and therefore B' is a 
maximum-weight base. So / -|- y is good. 

Sufficiency. Suppose that the greedy algorithm leads to an independent set of maximum 
weight for each weight function w : S ^ M+. We show that {S,I) is a matroid. 

Condition Equation 1.1 (i) is satisfied by assumption. To see condition Equation l.l(ii), 
let I,J£Z with |/| < \J\. Suppose that I + z ^ Z for each z £ J\I. Let k := |/|. Consider 
the following weight function w on S: 
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k + 2 if s e /, 
w{S) := I k + 1 Use J\I, 

ii s e S\{IU J). 
Now in the first iterations of the greedy algorithm we find the k elements in /. By 
assumption, at any further iteration, we cannot chose any element in J\I. Hence any 
further element chosen, has weight 0. So the greedy algorithm yields an independent set of 
weight k{k + 2). 

However, J has weight at least | J|(/c + 1) > (A; + l){k + 1) > k{k + 2). Hence the greedy 
algorithm does not give a maximum-weight independent set, contradicting our assumption. 

□ 



The rank function of a matroid Ai also determines Ai and we can easily see that a 
set ^ C S" is independent if and only if if {A) = \ A\. 

Theorem 6 (Theorem 39.8 [Sch03]). Let S be a set and ip : 2^ ^ Z+. Then Lp is the 
rank function of a matroid if and only if for all T,U S: 

(i) v{T)<ip{U)<\U\ifT<ZU, 

(1.4) 

(u) ip{TnU) + ip{TUU) <ip{T) + ip{U). 

Proof. Necessity. Let the rank function of a matroid {S,I). Choose T,U C S. 
Clearly Equation 1.4(i) holds. To see (ii), let I be an inclusionwise maximal set in X with 
I QTnU and let J be an inclusionwise maximal set in I with I Q J QTUU. Since (S, Z) 
is a matroid, we know that (p{T n f/) = |/| and (p{T UU) = \J\. Then 

ip{T)+ip{u) > |JnT|+| Jn[7| = |Jn(rn;7)|+| Jn(ruc/)| > \i\ + \J\ = ip{Tnu)+ip{Tuu); 

that is, we have Equation 1.4(ii). 

Sufficiency Let 2 be the collection of subsets / of 5 with (/?(/) = |/|. We show that 
{S,I) is a matroid, with rank function if. 

Trivially, G T. Moreover, if / G X and J G I, then 

^(J)<9'(/)-</^(/\J)>|/|-|/\J| = |J|. 
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So J e J. 

In order to check Equation 1.2, let I,J £ I with |/\ J| = 1 and |J\/| = 2. Let 
J\I = {zi,Z2}- li I + zi, I + Z2 ^ I, we have (/?(/ + zi) = ip{I + Z2) = \I\- Then by 
Equation 1.4(ii), 

ip{J) < ip{I + Zl + Z2) < ip{I + zi) + ip{I + Z2) - ip{I) = |/| < I J|, 

contradicting the fact that J 2. 

So {S,2) is a matroid. Its rank function is if, since ^piU) = max{ \I\ \ I <^ U, I £ 1} 
for each U £ S. Here > follows from Equation 1.4(ii), since if / C [/ and I £ I, then 
(/?([/) > (p{I) = 1-^1- Equality can be shown by induction on \ U\, the case C/ = being trivial. 
If [/ 7^ 0, choose y £ U. By induction, there is an / C [/ — y with I £ Z and |/| = Lp(U — y). 
If ^{U) = — y) we are done, so assume ip{U) > ip(U — y). Then I + y G X, since 
'^{I + y) > f{I) + <f{U)-^{U-y) > + Moreover, ip{U) < f{U-y) + ^{{y}) < |/| + 1. 
This proves equality for U. □ 

One of the most natural methods of presenting a matroid in a computational model 
is by independence oracle which given any A Q S asserts whether ^ G X or not. The 
independence oracle allows us to compute the rank of any A <^ S. Begin with / := and 
for each i £ A (any order) if / + z £ I then set I := I + i. Repeat and thus, rank(^) := |/|. 
For an interesting overview of the computational relation of various matroid axiomatization 
representations see [May08]. 

In Theorem 5 it was assumed that the weights are non-negative. But even with a weight 
function w : 5 — > M it can still be shown that for a matroid M the greedy algorithm finds a 
maximum- weight basis. Not surprising, one can replace the condition "as large as possible" 
with "as small as possible" and the greedy algorithm will find a minimum- weight basis. 
Assuming an independence oracle and Theorem 5 the corollary follows: 

Corollary 7 (Corollary 40.1a [Sch03]). A maximum- weight independent set in a ma- 
troid can be found in strongly polynomial time. 

Example 8 (Linear, Realizable, or Vector Matroids). Let A := {ai | 02 | ••• | £ 
jgmxn l)e an n X m matrix. Let S := {1, . . . ,n} and we say / C 5 is independent if the 



1.2. MATROID POLYTOPES & POLYMATROIDS 



7 



columns of A indexed by / are linearly independent. 

1 2 3 4 5 

( 1 1-1 2\ 
A = {1, 2}, {2, 3} are indepdendent while {1, 5} is not. 

\ 1 1 1 2 J 

Such a matriod Ma is called linear or realizable. If the entries of A are in a field F then 
we say Ma is realizable over F. Graphical matroids are linear and realizable over any field. 
This can be done by arbitrarily orienting the edges of G and take its incidence matrix. Not 
all matroids are linear, e.g. see Fg in the appendix of [0x192]. Some matroids can only be 
realized in certain fields. For example, the Fano matroid can only be realized in a field of 
characteristic two [0x192] . 

1.2. Matroid Polytopes & Poly matroids 

Let Ac [n] := {1, . . . ,n} and define the incidence vector of A as 

:= y^e^ 

where S are the standard unit vectors. For example if A = {1,3,4} C [5] then 
= (1,0,1,1,0)''' G M^. We also adopt the notation Inc(X) := {ba \ A £ X} where 
X C 2-^. Now we introduce the main object of this thesis. 

Definition 9. Let be a matroid with bases Bm- The matroid polytope of M is 

Vm ■= conv(eB | B G Bm) 

where conv(-) denotes the convex hull. 

Example 10. Recah U^'^ = { {1, 2}, {1, 3}, {1, 4}, {2, 3}, {2, 4}, {3, 4} }. Thus 

P[;2,4 =conv( (1,1, 0,0)"^, (1,0, 1,0)^,(1, 0,0, 1)^,(0, 1,1, 0)T, (0,1, 0,1)^,(0, 0,1, 1)T). 

The polytope P^2,4 G but we can easily see that all vertices lay in a hyperplane since 
they all sum to two. Vij2,i is in fact isomorphic to the cross polytope shown in Figure 1.2. 

The matroid polytope is different than the following well-known polytope. 
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(1,1,0,0) 




Figure 1.2. The cross polytope Vu2a. 

Definition 11. Let M = {S,I) be a matroid. The independence matroid polytope of 
M is 

Vm :=conv(e/ | /GT). 

We can see that V{M) C 'P-^(M) and Vm is a face of in the hyperplane X^j'Li Xi = 
rank(M). 

Polymatroids are closely related to matroid polytopes and independence matroid poly- 
topes. A function -0: 2^ — > M is suhmodular if ^(X n F) + ^{X \JY) < V'(X) + 
for all X, y C [n] and non- decreasing if V'(^) ^ V'l^) for all X C y C [n]. We say ip \s a 
polymatroid rank function if it is submodular, non-decreasing, and V'(O) = 0- For example, 
the rank function of a matroid is a polymatroid rank function. The polymatroid determined 
by a polymatroid rank function if) is the convex polyhedron in given by 

:= {x G M" I ^Xi< V'(^), C [n], x > o}. 

Independence matroid polytopes are a special class of polymatroids. Indeed, if (p is 
a rank function on some matroid A4, then Vj^ = "P,^ [Edm03]. Moreover, the matroid 
polytope Vm is the face of 'P^p lying in the hyperplane Yll=i ^« = '/'(N)- Matroid polytopes 
and polymatroids appear in combinatorial optimization [Sch03], algebraic combinatorics 
[FS05], and algebraic geometry [GGMS87]. Edmonds showed that a linear function can 
be optimized over a polymatroid via an extension of the greedy algorithm [Edm03]. 
This also provides a technique which easily lists off the vertices of a polymatroid V^. One 
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/ 

Figure 1.3. An example of a polymatroid 

important implication we will see is that if V' is an integral polymatroid rank function then 
will be integral, that is, all the vertices of are integral. 

First we recall some basic definitions and key theorems from linear programing [Sch86a] . 
If S" is a finite subset then M"^ is the vector space of vectors with coordinates indexed by S. 
That is, if X € M'^ and s ^ S then the sth coordinate of x is given by x{s) G M. li A C S 
then x(^) := Eae^a^l")- 

Let S, T be finite sets, c G , b G and A = {aijeR\ieT, j£S}hea matrix. 
A linear program is 

max{ c • X I Ax < b, X > 0, X G M"^ }. 

The dual linear program is 

min{ b • y I A^y < c, y > 0, y G }. 

A vector x G satisfying Ax < b is a feasible solution and a vector y G M.'^ satisfying 
< c is a feasible dual solution. 

Theorem 12 ([Sch03]). For any linear program maximization problem exactly one of 
the following is true; 

(i) There exists no feasible solution, 

(ii) For any q G M there is a feasible solution x such that c • x > a, 

(iii) There is an optimal (feasible) solution. 

A linear program and its dual are related and we state the famous strong duality theorem. 



1.2. MATROID POLYTOPES & POLYMATROIDS 



10 



Theorem 13 (Corollary 7.1g [Sch86a]). Let S,T be finite sets, AgR^ xET^, ce M'^, 
b G M'^. Then 

max{ c • X I < b, X > 0, X G R'^ } = min{ b • y | ^"^y < c, y > 0, y G M'^ } 

if both sets are nonempty. 

In order to justify the following approach to finding the vertices of a polymatroid we 
will use a fact from linear programming and polyhedral theory; every vertex of a polyhedra 
is the optimal solution to some linear program. Let <p he a polymatroid rank function on 
the finite set S and w : S" — > M. We start by ordering the elements of S" as si, . . . , s„ such 
that 'w{si) > w{s2) > ••• > w{sn) and define Ui := {si, . . . , Si} for z = 0, . . . , n. Next, 
define x G M'^ as 

(1.5) x{si) := ip{Ui) - ip{Ui-i) for i = 1, . . . , n. 

The following theorem will show that x maximizes w • x over V^p. Consider the linear 
program and its dual: 



(1.6) 



max{ w • X I X G Vip } 

^ y{T)ip{T) I y G , J] y(r) • > w 



mm < 



TCS TCS 



Also, define: 

y{Ui)= w{si) - w{si+i) (z = 1, . . . ,n - 1) 

(1-7) y{S)= w{sn), 

y(r) =0 {T^U, for each i). 

Theorem 14 (Theorem 44.3 [Sch03]). Let ip be a polymatroid rank function on S and 
let w : S ^ M*^. Then x and y are optimum solutions to Equation 1.6. 

Proof. We first show that x G Vip] that is, x < ip(T) for each T Q S. This is shown 
by induction on |r|, the case T = ^ being trivial. Let T 7^ and let k be the largest index 
with Sk G T. Then by induction. 



x(T-Sfc) < (^(T-Sfc). 
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Hence 

X(r) < v?(T - Sk) + X(xfe) = ip{T - Sk) + viUk) - viUk-i) < v{T) 

(the last inequality follows from the submodularity of (/?). So x G Vip. 

Also, y is feasible for Equation 1.6. Trivially, y > 0. Moreover, for any i we have by 
Equation 1.7: 

5^y(r)=^y(C/,)=w(.0. 

TBsi j>i 

So y is a feasible solution of Equation 1.6. 
Optimality of x and y follows from: 

n 

w • X = ^ w{s)x{s) = ^ ■w{si){ip{Ui) - ip{Ui-i)) 
ses 1=1 

n-l 

= ^'P{Ui){w{si) - w{si+i)) + ip{S)w{sn) = ^y{T)ip{T). 

i=l TC 

The third equality follows from a straightforward reordering of the terms, using that (/7(0) = 
0. □ 

Thus, Equation 1.5 gives a vertex of Vip. If we vary the weights we will get all the vertices 
of Vip and it suffices to choose weights that permute the elements of S. Equation 1.5 and 
Theorem 14 imply the following important corollary. 

Corollary 15 (Corollary 44.3d [Sch03]). Let (p be an integral polymatroid rank func- 
tion. Then the vertices ofVip are integral. 

Given a matroid M of rank r it is easy to see that Vm is contained in the (n — l)-simplex 

A = { (xi, . . . , x„) G M" I > 0, xi H \-Xn = r.} 

Moreover, the following theorem of Gel'fand, Goresky, MacPherson, and Serganova gives 
an equivalent geometric axiomatizatio of matroids: 

Theorem 16 ([GGMS87]). Let B be a family of r-element subsets of[n\. The following 
are equivalent: 

(i) B are the bases for a matroid, 

(ii) Every edge of conv^es \ B £ B) is parallel to an edge of the simplex A. 
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The following is a more useful characterization of adjacency of vertices of "Px, as well 
as a characterization of adjacency in the independence matroid polytope 

Lemma 17 (See Theorem 4.1 in [GGMS87], Theorem 5.1 and Corollary 5.5 in [Top84]). 
Let M he a matroid. 

(i) Two vertices e Bi and e B2 are adjacent in Vm and only if es^ ~ = — 
for some 

(ii) If two vertices e/^ and ej^ are adjacent in Vj^ then e/^ — e/j S {e^ — e^, e^, — } 
for some i,j. Moreover if \r is a vertex of then all adjacent vertices of \r can 
be computed in polynomial time in n, even if the matroid M is only presented by 
an evaluation oracle of its rank function (p. 

We prove one direction of (i), that every edge of Vm is parallel to — ej for some 
i,j G N. The following can be found in the proof of Proposition 2.2 in [Onn03]. 

Proof. Consider any pair A, B £ Bm of bases such that [6^1,65] is an edge (that is, 
a 1-face) of Vmi ^"^^ l^t a G be a linear functional maximized over Vm uniquely on 
that edge, li A\B = {i} is a singleton, then B\A = {j} is a singleton as well, in which 
case — = — and we are done. Suppose then, indirectly, that it is not, and 
pick and element i in the symmetric difference {A\B) {B \ A) of A and B of minimum 
value Oj. Without loss of generality assume i G A\B. Then there \s a j £ B\A such that 
C := A - i + j \s a, basis of M. Since \[A\B) U [B\A)\ > 2, C is neither A nor B. By 
the choice of i, this basis satisfies a • e^r = a • e^i — + > a • e^, and hence ec is also a 
maximizer of a over Vm ^-iid so lies in the 1-face [e/i,e5]. But no {0, l}-vector is a convex 
combination of other, yielding a contradiction. □ 

A similar property holds for integral polymatroids which we also state but we first recall 
some needed definitions from [Top84]. Let v,w G M" and define A(v,w) := { i G [n] | u j 7^ 
Wi} and cl(v) := { G | G C [n], Yli^c'^i ~ V'(G)}. Let F = {/i, • • • , /|_f|} be an ordered 
subset of [n] and Fj := {/i, . . . , /j}. If ■0 is a polymatroid rank function then we construct 
V G M" where Vi = ipiF'i) ~ V'(-^i-i) where vj = when j ^ F and one says F generates v. 

Lemma 18 (See Theorem 4.1 and Section 2 in [Top84].). Let ip be a polymatroid rank 
function. If v and w are vertices of the polymatroid V{ip) then either 
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(i) |A(v,w)| = 1 or 

(ii) cl(v) = cl(w) and A(v,w) = {c,d} for some c,d£ [n] where there exists some 
ordered set F = {/i, . . . ,/|f|} which generates v with fk+i = d and fk = c for 
some integer k, 1 < A; < |i^| — 1; moreover the ordered set 

F ■= {fi,---, /fc-i, /fc+i, fk, fk+2, f\F\} generates w. 

In order to describe the matroid polytope Vm by a system of linear inequalities we 
require a few more definitions. A circuit C C 5 of = {S,I) is an inclusionwise minimal 
dependent set. That is, a dependent set C is a circuit of M if there does not exist a 
dependent set A such that A C C Q S. A flat F of is a subset of S such that 
ip{F + s) = ip{F) + 1 for all s £ S. The span of A C 5 is the smallest flat F such that 
A <^ F. A matroid is called connected if for every x,y G S there exists a circuit C such that 
x,y £ C. In fact the circuits of a matroid M define an equivalence class on the ground set 
S where x, y are equivalent if there exists a circuit C such that x,y G C. The connected 
components of M are defined as the number of such equivalence classes. 

Lemma 19 ([FS05]). The matroid polytope Vm equals the following subset of the simplex 

A: 

'Pm = \ (2^1 > • • • ) Xn) G A I Xi < rank(F) for all fiats F <^ S > 
{ ieF ) 

Proof. Consider any facet of the polytope Vm and let Yll=i '^i^i < ^ be an inequality 
defining this facet. The normal vector (oi, 02, ... , a„) is perpendicular to the edges of that 
facet. But each edge of tht facet is parallel to some difference of unit vectors — Cj by 
Lemma 17. Hence the only constraints on the coordinates of the normal vector are of the 
form ai = aj. Using the equation X^ILi -^^ ~ ^^"^ scaling the right hand side b, we can 
therefore assume that (oi, 02, • • . , a^) is a vector in {0,1}". Hence the polytope Vm is 
characterized by the inequalities of the form ^^^q Xi < be for some G Q [n]. The right 
hand side be equals 

bo = max{ |(t| n G I (T basis of M} = rank(G). 

The second equality holds because every independent subset of G can be completed to a 
basis a. Let F be the flat spanned by G. Then G Q F and rank(G) = rank(F), and hence 
the inequality "^^^pXi < rank(F) implies the inequality "^^^QXi < rank(G). □ 
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Lemma 20 ([FS05]). The dimension of the matroid polytope Vm equals n — ^connected 
components of M. 

Proof. Two elements i and j are equivalent if and only if there exist bases a and r 
with i £ a and r = a + i—j. The linear space parallel to the affine span of Vm is spanned by 
the vectors — ej arising in this manner. The dimension of this space equals n — c(M). □ 

1.2.1. Matroid Optimization. Matroids find application throughout combinatorics. 
The matroid intersection problem is defined as finding the largest size common independent 
set of two matroids. We have seen that the greedy algorithm fully characterizes matroids 
and Edmonds proved that matroid intersection, can be solved efficiently. This implies that 
bipartite matchings, common transversals, and tree packing and covering can be computed 
efficiently too [Sch03] . Interestingly, it has been shown that finding matroid intersection of 
three or more matroids is NP-complete. 

Matroids are particularly important in optimization. The greedy algorithm requires 
weightings on the ground set S of the matroid and it has many applications to optimization. 
But in real life, multiple weightings can be applied to the ground set along with a balancing 
function, or concrete notion of optimality over multiple weightings. Consider d weightings 
wi,...,W(i G M" on the ground set [n]. That is, every Wj assigns a real- value to each 
element of [n]. We let W E M'^^" be the matrix with rows wi, . . . , w^. We also define 
WVm := {WeB \ B € Bm} #WVm to mean the cardinality of WVm, WVm ■= 

conv(WVM) Slid vert('P) = {vertices of V} for a polytope V. Intuitively one thinks of the 
rows of W as the set of criteria that (possibly confiicting) parties may bring to a discussion. 
For example, think of two groups (loggers vs. environmentalists) trying to decide how to 
choose an optimal spanning tree (a network of logging roads) of a graph. Each group will 
have different weights. 

Hence in Chapter 4 of this dissertation, we will consider four generalizations of the 
classical single-criterion (linear-objective) matroid optimization problem (i.e. d = 1 which 
is well-known to be solvable via the greedy method). 

Non-Linear Matroid Optimization: Given a matroid M on [n] with bases i3x, 
W E M*^^", and a function / : M'^ ^ M, find a base B e Bm such that fiWes) = 
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mm{ f{WeB') \ B' e B). 

Two important special cases of our investigations are the following problem. 

Convex Matroid Optimization: Given a matroid Ai on [n] with bases B, W ^ W^^"^ 
where d is fixed, and a convex function / : M*^ ^ M, find a base B ^ B such that 
fiyVes) = max/(l/V^e^/ | B' G B). Similarly we consider the minimization problem 
fiWee) = min /(l^e^/ \ B' e B) 

Min-Max Multi-criteria Matroid Optimization: Given a matroid M. on [n] with 
bases Bm, W G M°'^", find a base B G Bm such that 

max {{WeB)i) = min( max {{WeBi)i) \ B' G ). 

i=l...d i=l...d 

Also, we investigate the following problem. 

Pareto Multi-criteria Matroid Optimization: Given a matroid M on n-elements 
with bases Bm, W G M'^^'^, find a base B G Bm such that 
S = argminp„^gto( (VFeij/) | G )• 

In the previous statement, minpa^eto is understood in the sense of Pareto optimality 
for problems with multiple objective functions, namely, we adopt the convention that for 
vectors a, b G M"', we have a < b if and only if aj < 6j for all entries of the vectors. Further, 
we say that a < b if a < b and a 7^ b. We note that an optimum of the minmax problem 
will be a Pareto optimum. This is easy to see since if a < b then maxi(aj) < maxj(6j). The 
Pareto multi-criteria matroid optimization problem has been studied by several authors 
before. For example, Ehrgott [Ehr96] investigated two optimization problems for matroids 
with multiple objective functions, and he pioneered a study of Pareto bases via the base- 
exchange property of matroids. See [EGOO] for a detailed introduction to multicriteria 
combinatorial optimization. 

The matroid optimization problems we consider here have wide applicability. For exam- 
ple, in [BLMA+08] the authors consider the "minimum aberration model fitting" problem 
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in statistics, which can be reduced to a non-hnear matroid optimization problem. Multi- 
criteria problems concerning minimum spanning trees of graphs are common in applications 
(see [EGOO] and [KC02] and references there in). 

The Min-Max optimization problem includes the NP-complete partition problem [G J 79], 
certain multi- processor scheduling problems [GLLR79], and specific worst-case stochastic 
optimization problems [War85]. 

Example 21 (Multi-criteria Spanning- Tree Optimization). Consider a graph G with 
several linear criteria on the edges, 

1. the first row of W encodes the fixed installation cost of each edge of G; 

2. the second row of W encodes the monthly operating cost of each edge of G; 

3. assuming that the edge j fails independently with probability 1 — pj, then by 
having the logpj as the third row of W (scaled and rounded suitably), X^jg^'logpj 
captures the reliability of the spanning tree T of G. 

It can be difficult for a decision maker to balance these three competing objectives. There 
are many issues to consider such as the time horizon, repairability, fault tolerance, etc. 
These issues can be built into a concrete function /, for example a weighted norm, or can 
be thought of as determining a black-box /. 

Unfortunately, although useful, the problems we are considering are also very difficult 
in general. For example, multi-criteria matroid optimization is generally NP-hard even in 
the case when d = 2 for uniform matroids. It is also NP-hard for the case of spanning trees. 
Nevertheless, recently many theoretical and complexity properties about these problems 
have been explored, yielding good complexity bounds and algorithms under nice assump- 
tions about W, f and d. For instance, it has been shown that although convex matroid 
optimization is NP-complete in general, it is polynomial-time solvable under certain restric- 
tions on W and with fixed d. We refer the reader to the recent series of papers on nonlinear 
matroid optimization [BLOW08, BOOS, Onn03, BLMA+08]. This dissertation explored 
the experimental performance of their algorithms and proposed some new heuristics and 
algorithms which performed well in practice. We will see this in Chapter 4. 
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CHAPTER 2 

Volume and Ehrhart Polynomial Computation 

To state our main results recall that given an integer A; > and a polytope V C M" 
we define kV := { kcx | a G "P } and the function i('P, k) := ij^{kV H Z"), where we define 
i('P, 0) := 1. It is well known that for integral polytopes, as in the case of matroid polytopes, 
i('P, k) is a polynomial, called the Ehrhart polynomial of V. Moreover the leading coefficient 
of the Ehrhart polynomial is the normalized volume of "P, where a unit is the volume of the 
fundamental domain of the affine lattice spanned by V [Sta96]. Our first theorem states: 

Theorem 22. Let r he a fixed integer. Then there exist algorithms whose input data 
consists of a number n and an evaluation oracle for 

(a) a rank function (p of a matroid M on n elements satisfying '^{A) < r for all A, or 

(b) an integral polymatroid rank function ip satisfying tp{A) < r for all A, 

that compute in time polynomial in n the Ehrhart polynomial (in particular, the volume) 
of the matroid polytope V{M), the independence matroid polytope V^{M), and the polyma- 
troid V (ip) , respectively. 

The computation of volumes is one of the most fundamental geometric operations and 
it has been investigated by several authors from the algorithmic point of view. While there 
are a few cases for which the volume can be computed efficiently (e.g., for convex polytopes 
in fixed dimension), it has been proved that computing the volume of polytopes of varying 
dimension is #P-hard [DF88, BW91, Kha93, LawDla]. Moreover it was proved that even 
approximating the volume is hard [Ele86]. Clearly, computing Ehrhart polynomials is a 
harder problem still. To our knowledge two previously known families of varying-dimension 
polytopes for which there is efficient computation of the volume are simplices or simple poly- 
topes for which the number of vertices is polynomially bounded (this follows from Lawrence's 
volume formula [Law91a]). Already for simplices it is at least NP-hard to compute the whole 
list of coefficients of the Ehrhart polynomial, while recently [Bar06] presented a polynomial 
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time algorithm to compute any fixed number of the highest coefficients of the Ehrhart poly- 
nomial of a simplex of varying dimension. Theorem 22 provides another interesting family 
of varying dimension with volume and Ehrhart polynomial that can be computed efficiently. 
The proof of Theorem 22, presented in Section 2.1, relies on the geometry of tangent cones 
at vertices of our polytopes as well as a new refined analysis of the evaluation of Todd 
polynomials in the context of the computational theory of rational generating functions 
developed by [Bar94, BP99, BW03, Bar06, DLHTY04, DLHH+04, Woo04, VW08]. A nice 
introduction to these topics can be found at [BR07] 



2.1. Computing the Ehrhart Polynomials 

2.2. Preliminaries on Rational Generating Functions. 

Generating functions are crucial to proving our main results. For a good reference for the 
basic concepts used here see [BR07]. Let V C M" be a rational polyhedron. The multivariate 
generating function ofV is defined as the formal Laurent series in ^[[zi, . . . , Zn, z^^ , . . . , z~^] 

5p(z)= Yl 

where we use the multi-exponent notation = ]Xi=i -^f ' • If 'P is bounded, g-p is a Laurent 
polynomial, which we consider as a rational function g-p. liV is not bounded but is pointed 
(i.e., V does not contain a straight line), there is a non-empty open subset U C C" such that 
the series converges absolutely and uniformly on every compact subset of C/ to a rational 
function gp (see [BP99] and references therein). If V contains a straight line, we set gp = 0. 
The rational function gp G Q(zi, . . . , z„) defined in this way is called the multivariate 
rational generating function of P Pi Z". Barvinok [Bar94] proved that in polynomial time, 
when the dimension of a polyhedron is fixed, gp can be represented as a short sum of 
rational functions 

where £i e {—1, 1}. 

Our first contribution is to show that in the case of matroid polytopes of fixed rank, 
this still holds even when their dimension grows. Let v be a vertex of V. Define the tangent 
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cone or supporting cone of v to be 

C-p(v) := { V + w I V + ew e V for some e > } . 

We rely on the following result, which connects the rational generating function of a rational 
polyhedron to those of the tangent cones of its vertices. This result was independently 
discovered by Brion [Bri88] and Lawrence [Law91b]. A proof can also be found in [BP99] 
and [BHS09]. 

Lemma 23 (Brion-Lawrence's Theorem) . Let V he a rational polyhedron and V{V) he 
the set of vertices ofV. Then, 

where Cp(v) is the tangent cone ofy^. 

Thus, we can write the multivariate generating function of V by writing all multivariate 
generating functions of all the tangent cones of the vertices of V. Moreover, the map 
assigning to a rational polyhedron V its multivariate rational generating function g-pi^) is 
a valuation, i.e., a finitely additive measure, so it satisfies the equation 

for arbitrary rational polytopes Vi and V2, the so-called inclusion-exclusion principle. This 
allows to break a polyhedron V into pieces Vi and and to compute the multivariate 
rational generating functions for the pieces (and their intersection) separately in order to 
get the generating function g-p. More generally, let us denote by [V] the indicator function 
of V, i.e., the function 



[P]:M"^M, [P](x) 



f 

1 if X e 

otherwise. 



Let X]jG/^«[^«] = be an arbitrary linear identity of indicator functions of rational poly- 
hedra (with rational coefficients e^); the valuation property now implies that it carries over 
to a linear identity Yli^i dViiT') = of rational generating functions. 
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Now let C be one of the tangent cones of V, and let T be a triangulation of C, given by 
its simplicial cones of maximal dimension. Let T denote the set of all (lower-dimensional) 
intersecting proper faces of the cones Ci £ T. Then we can assign an integer coefficient Ei 
to every cone Ci E T, such that the following identity holds: 

[C]= Et^^l+E^^P^l- 
This identity immediately carries over to an identity of multivariate rational generating 
functions, 

(2.1) 5c(z) = Yl 9cA^) + E ^*5c,(z). 

c.e^ c.er 

Hence, the problem of computing rational generating functions of a polyhedron is reduced 
to the case of simplicial cones. 

In the following (section 2.3, section 2.4), we study the tangent cones of the matroid poly- 
topes and polymatroids and introduce algorithms that construct triangulations for them. 
Then, in section 2.5, we construct a short multivariate generating function using an efficient 
variant of identity (2.1) and Brion-Lawrence's Theorem. Finally, in section 2.6, we compute 
the Ehrhart polynomial. 

2.3. On the Tangent Cones of Matroid Polytopes 

Our goal is to compute the multivariate generating function of matroid polytopes and 
independence matroid polytopes with fixed rank (later, in section 2.4, we will deal with the 
case of polymatroids), and to do this we will use a crucial property of adjacent vertices. To 
illustrate our techniques we will use a running example throughout this section. 

Example 24 (Matroid on K4). Let K4 be the complete graph on four vertices. Label 
the (2) = 6 edges with {!,..., 6} as in Figure 2.1. Every graph induces a matroid on 
its edges where the bases are all spanning trees (spanning forests for disconnected graphs) 
[Wel76]. Let M(K4) be the matroid on the elements {1, . . . , 6} with bases as all spanning 
trees of K4. The rank of M(i^4) is the size of any spanning tree of -R'4, thus the rank of 
M(i^4) is 3. The 16 bases M^K^) are: {3,5,6}, {3,4,6}, {3,4,5}, {2,5,6}, {2,4,6}, 
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{2,4,5}, {2,3,5}, {2,3,4}, {1,5,6}, {1,4,6}, {1,4,5}, {1,3,6}, {1,3,4}, {1,2,6}, {1,2,5}, 
{1,2,3}. 

Let M be a matroid on n elements with fixed rank r. Then the number of vertices of 
V{M) is a polynomial in n of degree r. We can see this since the number of vertices is equal 
to the number of bases of M, and the number of bases is bounded by , a polynomial in 
n of degree r. Clearly the number of vertices of V^{M) is also polynomial in n. It is also 
clear that, when the rank r is fixed, all vertices of either polytope can be enumerated in 
polynomial time in n, even when the matroid is only presented by an evaluation oracle for 
its rank function ip. 

Throughout this section we shall discuss polyhedral cones C with extremal rays {ri, . . . , r^} 
such that 

Tfc e := {ej - ej, e^, -e^- | i £ [n], j e A} for = 1, . . . , ^ 

for some A CI [n]. We will refer to Ra as the elementary set of A. Note that by Lemma 17 the 
rays of a tangent cone at a vertex e^^ (corresponding to a set A C [n]) of a matroid polytope 
or an independence matroid polytope form an elementary set of A. Due to convexity and the 
assumption that are extremal, for each i £ [n] and j £ A at most two of the three vectors 
ej — Cj, ej, — ej are extremal rays of C. This implies by construction, that considering all 
pairs Cj — ej and or —ej, the number of generators of C is bounded by 

(2.2) n\A\+n + \A\. 

Recall a cone is simple if it is generated by linearly independent vectors and it is uni- 
modular if its fundamental parallelepiped contains only one lattice point [G097]. A tri- 
angulation of C is unimodular if it is a polyhedral subdivision such that each sub-cone is 
unimodular. 

Example 25 (Matroid on K^). The vertices e{2,3,5}, e{2,3,4}, eji^gg}, e{i^34}, e{i^2,6} 
and e|i 2,5} are all adjacent to the vertex ej^ 2,3}) see Figure 2.1. Moreover, the tangent 
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Figure 2.1. {2,3,5}, {2,3,4}, {1,3,6}, {1,3,4}, {1,2,6} and {1,2,5} are 
spanning trees of K4, that differ from {1,2,3} by adding one edge and re- 
moving one edge. 

cone C-p(jvf (K4))(e{i^2,3}) is generated by the differences of these vertices with e|i 2,3}: 

Cp{M(K4))(6{l,2,3}) = ^{1,2,3} + cone {e|2,3,5} - e{i_2,3};e{2,3,4} " ^{1,2,3}! 

^{1,3,4} - e{i,2,3}ie{i_2,6} - 2,3}>e{i,2,5} -^{1,2,3}} 

Lemma 26. Let C Q M*^ be a cone generated by p extremal rays {ri, . . . , Vp} C i?^ where 
Ra is an elementary set of some A CI [n] . Every triangulation of C is unimodular. 

Proof. Without loss of generahty, we can assume {ri, . . . ,ti\ are generators of the 
form ej — and {r^+i, . . . , Yp} are generators of the form ej or — for the cone C. 

It is easy to see that the matrix Tq := [ri, . . . , r;] is totahy unimodular. Let Gq be a di- 
rected graph with vertex set [n] and an edge from vertex i to j if = — Gj is an extremal 
ray of C. We can see that Gq is a subgraph of the complete directed graph Kn with two arcs 
between each pair of vertices; one for each direction. Since Tq := [ri, . . . , r;] is the incidence 
matrix of the graph Gq, it is totally unimodular [Sch86b, Ch. 19, Ex. 2], i.e., every subde- 
terminant is 0, 1 or —1 [Sch86b, Ch. 19, Thm. 9]. Therefore Tq := [ri, . . . , r;, r^+i, . . . , rp] 
is totally unimodular since augmenting Tq by a vector or — preserves this subdetermi- 
nant property: for any submatrix containing part of a vector or —ej perform the cofactor 
expansion down the vector or — e^- when calculating the determinant. 
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Since Tq is totally unimodular, each basis of Tc generates the entire integer lattice 
Z" nlin(C) and hence every simplicial cone of a triangulation has normalized volume 1. □ 

Lemma 27. Let C Q be a cone generated by I extremal rays {ri, . . . , r^} C Rj^ where 
Ra is an elementary set of some A CI [n], where dim(C) < n. The extremal rays {ri, . . . , r^} 
can be augmented by a vector r such that dim(cone{ri, . . . , r/, f}) = dim(C) + 1, the vectors 
ri, . . . , r;, f are all extremal, and f G Ra- 

Proof. It follows from convexity that at most two of ej — ej, ej or —ej are extremal 
generators of C for i S [n] and j E A. There are at least n possible extremal ray generators, 
considering two of — e,,-, ej or — e,,- for each i G [n] and j € A. Moreover, all these pairs 
span M". Thus by the basis augmentation theorem of linear algebra, there exists a vector f 
such that dim(cone{ri, . . . , r^, f}) = dim(C) + 1 and ri, . . . , r^, f are all extremal. □ 

Lemma 28. Let r be a fixed integer, n be an integer, A C [n] with \ A\ < r and let C C M" 
be a cone generated by I extremal rays {ri, . . . , r^} C Ra where Ra is an elementary set of 
A. Then any triangulation of conv({0, ri, . . . , r;}) has at most a polynomial in n number 
of top- dimensional simplices. 

Proof. Assume dim(C) = n. Later, we will show how to remove this restriction. We 
can see that conv{0, ri, . . . , r;} C [—1, 1]^ x A[„]\^ where 

[-1,1]^ := {x e I |xj| < 1 jeyljcM^ 

A[„,]\^ := conv({ei | i G [n] \ A} U {0}) C mM\^. 
The volume of a d-simplex convjvo, . . . , v^^} is [G097] 

detf "° \ . 

\l ... 1 j 

Thus the (n — |A|)-volume of A[„]\^^ is ^.^_|^|^| and the |^|-volume of [—1, 1]"^ is 2l^L There- 
fore 

vol ([-1, 1]^ X A[„]\^) = 2l^l (^^3^ = l,2\Mn - 1) • • • (n - + 1). 

It is also a fact that any integral n-simplex has n- volume bounded below by ^, using the 
simplex volume equation (2.3). Therefore any triangulation of conv{0, ri, . . . , r/} has at 



(2.3) 



1 



2.3. ON THE TANGENT CONES OF MATROID POLYTOPES 



24 



most 

2l^ln(n-l)---(n-|A|+l) < 2^ n(n — 1) • • • (n — r + 1) 

full-dimensional simplices, a polynomial function in n of degree r. 

Let dc := n — dim(C). If dim(C) < n, then by Lemma 27, {ri, . . . ,r/} can be aug- 
mented with vectors {f i, . . . , fd^} where ffc G Ra for A above, such that dim(cone{ri, . . . , r;, 
fi, . . . , Ydc)) = ''^ ^-iid {ri, . . . , r^, f 1, . . . , f^^} are extremal. Moreover, 

dim(conv{0, ri, . . . , r^}) < dim(conv{0, ri, . . . , r^, f i}) < • • • 
< dim(conv{0,ri, . . . ,r;,fi, . . . ,fd^_i}) < dim(conv{0, ri, . . . ,ri,fi, . . . jf^^}), 

that is, ffc i afr{0, ri, . . . , r/, f i, . . . , ffc_i} for I <k <dc. 




Figure 2.2. If is not contained in the affine span of 

{0, ri, . . . , Tp, f 1, . . . , ffc-i} then every full-dimensional simplex must 
contain f^. 



Since ffc ^ aff{0, ri, . . . , r^, f i, . . . , f fc-i}, any full- dimensional simplex in a triangulation 
of convjO, ri, . . . , r/, f 1, . . . , r^} must contain f^, see Figure 2.2. If not, then there exists a 
top-dimensional simplex using the points {0, ri, . . . , r;, ri, . . . , ffe_i}, but we know all these 
points lie in a subspace of one less dimension, a contradiction. Therefore, a bound on the 
number of simplices in a triangulation of convjO, ri, . . . , r^, fi, . . . , r^} is a bound on that 
of conv{0,ri, . . . ,r;,fi, . . . ,ffc-i}- 

Thus, if dim(C) < n we can augment C by vectors fi, . . . ,fdg so that the cone C := 
conejri, . . . , r^, fi, . . . , ff^^} is of dimension n and r^, f^. G Ra for A above. We proved 
any triangulation of convjO, ri, . . . , r^, f i, . . . , frf^j} has at most polynomially many full- 
dimensional n-simplices, which implies that any triangulation of conv{0, ri, . . . , Yp} has at 
most polynomially many top-dimensional simplices due to the construction of the generators 
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We have shown that for a cone C generated by an elementary set of extremal rays 
{ri, . . . ,r/} C for some A C [n], any triangulation of conv{0,ri, . . . ,r;} has at most 
polynomially many simplices. What we need next is an efficient method to compute some 
triangulation of convjO, ri, . . . , r;}. We will show that the placing triangulation is a suitable 
candidate. 

Let V C M" be a polytope of dimension n and A be a facet of V and v G M" . There 
exists a unique hyperplane H containing A and V is contained in one of the closed sides of 
H, call it H^. If V is contained in the interior of H~ , the other closed halfspace defined 
by H, then A is visible from v (see chapter 14.2 in [G097]). The well-known placing 
triangulation is given by an algorithm where a point is iteratively added to an intermediate 
triangulation by determining which facets are visible to the new point [G097, DLRS09]. 
We recall now how to determine if a facet is visible to a vertex in polynomial time. 




Lemma 29. Let V QW^ be a polytope given by t vertices {v^, . . . , v*} C M" and A C P 
be a facet of V given by q vertices {^r^, ■ ■ ■ , v*^} ^ {v^, . . . , v*}. If v £ R" where v ^ "P 
then deciding if A is visible to v can be done in polynomial time in the input {v^, . . . , v'^}, 
{v-*^, . . . , V*} and v. 
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Proof. Let z := - Yli=i that z S relint(A). We consider the Hnear program: 




t 



y 



n+t+l 



X = ^vVi,y 



>o, 5^y. = 1 



(2.4) 



1=1 



< A < 1, Av + (1 - A)z 



= X 



If (2.4) has a solution then there exists a point x G "P between the facet A and v, hence 
A is not visible from v. If (2.4) does not have a solution, then there are no points of V 
between v and A, hence A is visible from v (see Lemma 4.2.1 in [DLRS09]). It is well 
known that a strict inequality, such as the one in (2.4), can be handled by an equivalent 
linear program which has only one additional variable. Determining if (2.4) has a solution 
can be done in polynomial time in the input [Sch86b]. □ 



The placing triangulation is obtained by incrementally adding one point at a time, 
connecting the new point to the current triangulation. More precisely: 
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Algorithm 30 (The Placing Triangulation [G097, DLRS09]). 
Input: A set of ordered points {vi, . . . , vj} G M". 
Output: A triangulation T of {vi, . . . , vt} 
1: T:={{vi}}. 

2: for each Vj G {v2, . . . , vj} do 

3: Let B £T. 

4: := {vi, . . . , Vi_i} 

5: if Vi ^ aff(Pi) then 

6: T := 

7: for each D G T do 

8: T ■.= T' \j{D\j{wi}}. 

9: else 

10: for each B & T and each {\B\ — l)-subset C of i? do 

11: Create and solve the linear program (2.4) with (Pj, C, Vj) to decide 

visibility of C to Vj. 
12: if C is visible to Vj then 

13: T' :=T'u{C7U{vi}} 

14: T := T 
15: return T 

Indeed, Algorithm 30 returns a triangulation [G097]. We will show that for certain 
input (a point set corresponding to a vertex cone of a matroid polytope or independence 
matroid polytope), it runs in polynomial time. We remark that there are exponentially, 
in n, many lower dimensional simplices in any given triangulation. But, it is important to 
note that only the highest dimensional simplices are listed in an intermediate triangulation 
(and thus the final triangulation) in the placing triangulation algorithm. 



Theorem 31. Let r he a fixed integer, n he an integer, A C [n] with \A\ < r, and let 
{ri, . . . , r^} C Rj^. Then the placing triangulation (Algorithm 30) with input {0, ri, . . . , r^} 
runs in polynomial time. 
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Proof. By Equation (2.2) there is only a polynomial, in n, number of extremal rays 
{ri, . . . , Ti}. Thus, the for statement on line 2 repeats a polynomial number of times. Step 
5 can be done in polynomial time by solving the linear equation [ri, . . . , ri_i]x = v^. 

The for statement on line 7 repeats for every simplex D in the triangulation T, and the 
number of simplices in T is bounded by the number of simplices in the final triangulation. 
By Lemma 28 any triangulation of extremal cone generators in with the origin will use at 
most polynomially many top-dimensional simplices. Hence the number of top-dimensional 
simplices of any partial triangulation T will be polynomially bounded since it is a subset of 
the final triangulation. 

The for statement on line 10 repeats for every simplex B and every {\B\ — l)-simplex of 
B. As before, the number of simplices B is polynomially bounded, and there are at most n 
{\B\ — l)-simplices of B. Thus the for statement will repeat a polynomial number of times. 

Finally, by Lemma 29, determining if C is visible to Vj can be done in polynomial time. 
Therefore Algorithm 30 runs in a polynomial time. □ 



Corollary 32. Let r be a fixed integer, n be an integer, A C [n] with \A\ < r, and let 
C C M" be a cone generated by extremal rays {ri, . . . ,ri} C Rj^. A triangulation of C can 
be computed in polynomial time in the input of the extremal ray generators {ri, . . . ,r/}. 



Proof. Let Vc ■= convjO, ri, . . . , r^}. We give an algorithm which produces a trian- 
gulation of Vc '■= convjO, ri, . . . , rt} such that each full-dimensional simplex has as a 
vertex. Such a triangulation would extend to a triangulation of the cone C. This can be 
accomplished by applying two placing triangulations: one to triangulate the boundary of 
Vc not incident to 0, and another to attach the triangulated boundary faces to 0. The 
algorithm goes as follows: 

1) Triangulate Vc using the placing triangulation algorithm. Call it T'. 

2) Triangulate Vc using the boundary faces of T' which do not contain v. 
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Algorithm 33 (Triangulation joining to boundary faces) . 
Input: A triangulation T' of "Pc, given by its vertices. 
Output: A triangulation T of Vc such that every highest dimension 
simplex of T is incident to 0. 
1: T := 

2: for each C where C is a (|-B| — l)-simplex oi B ^T' do 

3: if C is not a (|^| - l)-simplex oi AeT' where B then 

4: r:=ru{cu{o}} 

5: return T 



Figure 2.4. A triangulation T' of Vc can be used to extend to a triangu- 
lation T such that is incident to every highest dimensional simplex 

By Theorem 31, triangulating Vc using Algorithm 30 can be done in polynomial time. 
Algorithm 33 indeed produces a triangulation of Vc- It covers Vc since every extremal ray 
generator r^ of C is on some (dim(C) — l)-simplex. Moreover, T by construction has the 
property that the intersection of any two simplices of T is a simplex. Step 3 checks if C 
is on the boundary, since if C is on the boundary it will not be on the intersection of two 
higher-dimensional simplices. 

Step 2 repeats a polynomial number of times since any triangulation of Vc has at most 
a polynomial number of simplices, and each simplex B has at most n {\B\ — l)-simplices. 
Step 3 can be computed in polynomial time since again there are only polynomially many 
simplices B in the triangulation T' and at most n {\B\ — l)-simplices to check if they are 
equal to C. Hence, Algorithm 33 runs in polynomial time. □ 



(a) 




(b)"b 
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Example 34. The tangent cone at the vertex := ej^ 2,3} on the polytope V{M{K4)) 
can be triangulated as: 

I {^{2,3,5} - ^B, e{2,3,4} " ^B, ^{1,^,6} - ^B, ^{1,3,4} - 6^, e{i^2,6} - ^b} , 

{^{2,3,5} - ^B, e{i,3,6} - e^, e{i^3_4} - cb, e{i_2,6} - e^, e{i^2,5} - ^b] , 
{^{2,3,5} - es, e{2,3,4} - e^, e{i^3_4} - e^, e{i 2,6} - e^, e{i^2,5} - ^b] |- 

2.4. Polymatroids 

We will show that certain lemmas from Subsection 2.3 also hold for certain polymatroids. 
Recall that the rank of the matroid M is the size of any basis of M which equals (^([n]). Our 
lemmas from Subsection 2.3 rely on the fact that M has fixed rank, that is, for some r G Z, 
^ ^ 0, '^{A) < r for all A CI [n]. We will show that a similar condition on a polymatroid 
rank function is sufficient for the lemmas of Subsection 2.3 to hold. 

Lemma 35. Let ip: 2^"'^ — > N be an integral polymatroid rank function where ^{A) < r 
for all ACL [n], where r is a fixed integer. Then the number of vertices ofV^ip) is bounded 
by a polynomial in n of degree r. 

Proof. It is known that if ip is integral then all vertices of 'P(V') are integral [Wel76]. 
The number of vertices of V{'ip) can be bounded by the number of non-negative integral 
solutions to + • • • + < r, which has ("^'') solutions, a polynomial in n of degree r 
[Sta97]. □ 

Lemma 36. Let 2["1 — > N be an integral polymatroid rank function. If v is a vertex 
of V{'4)) then all adjacent vertices of v can be enumerated in polynomial time. Moreover 
if ip{A) < r for all A C [n], where r is a fixed integer, then the vertices of V{ij)) can be 
enumerated in polynomial time. 

Proof. If v is a vertex of V{ip) then generating and listing all adjacent vertices to v 
can be done in polynomial time by Lemma 17. If ip{A) < r for all A Q [n], where r is a 
fixed integer, then, by Lemma 35, there is a polynomial number of vertices for V^ip). We 
know that S M" is a vertex of any polymatroid. Therefore, beginning with 0, we can 
perform a breadth-first search, which is output-sensitive polynomial time, on the graph of 
V{ip), enumerating all vertices of V^ip). □ 
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What remains to be shown is that these polymatroids have cones Hke the ones in Sub- 
section 2.3. 

Lemma 37. Let ip be an integral polymatroid rank function and C the tangent cone of a 
vertex v of the polymatroid V{ip), translated to the origin. Then C is generated by extremal 
ray generators {ri, ■ ■ ■ ,ri} C Rsupp{v)7 where i?supp(v) ^-5 an elementary set o/supp(v). 

Proof. Let tp: 2["'l — > Z be a integral polymatroid rank function. Let v and w 
be adjacent vertices of the polymatroid V{'ip). Using Lemma 18, if |A(v,w)| = 1 then 
w — V = hei where h is some integer and is the standard ith elementary vector for some 
i G [n]. If /i < then certainly i E supp(v), else i G [n]. Thus w — v, a generator of C, is 
parallel to a vector in iisupp{v)- 

Let V and w be adjacent and satisfy (ii) in Lemma 18, where A(v,w) = {c, d}. Hence 
there exists an F = {/i, . . . , /|_p|} which generates v with f^^i = d and fk = c for some inte- 
ger A;, 1 < A; < |F| - 1; moreover the ordered set F := {fi, . . . , fk-i, fk+i, fk, fk+2, f\F\} 
generates w. First we note that ■ip{Fk-i) = tp{Fk-i) and ip{Fk+i) = ^p{Fk+l)■ By assump- 
tion, we know Vc / Wc, Vd / Wd and vi = wi for all / G [n] \ {c, d}. Thus 

{^r-w)c = Vc-Wc = tP{Fk+i) - i>{Fk) - (i>{Fk) - ij{Fk^i)) 

= i^{Fk+i) - i>{Fk) - 4^{Fk) + ^{Fk-i) 

and 

(v - w)d = Vd-Wd = i^{Fk) - i>{Fk-i) - (i>{Fk+i) - i^{Fk)) 

= i;{Fk) - HFk-i) - (^(Ffc+i) - V(Ffc)) 
= -ij{Fk+i) + iP{Fk) + iP{Fk) - i^{Fk-i). 

Therefore (v — w)c = — (v — w)^ and w — v is parallel to — ec. Moreover, c G supp(v) 
since w, v > by assumption that v, w G V{tp). Thus w — v, a generator of C, is parallel 
to a vector in -Rsupp(v)- D 
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2.5. The construction of a short multivariate rational generating function 

From the knowledge of triangulations of tangent cones of matroid polytopes, indepen- 
dence matroid polytopes, and polymatroids we will now recover short multivariate gener- 
ating functions. 

Remark 38. Notice that formula (2.1) is of exponential size, even when the triangu- 
lation T only has polynomially many simplicial cones of maximal dimension. The reason 
is that, when the dimension n is allowed to vary, there are exponentially many intersect- 
ing proper faces in the set T. Therefore, we cannot use (2.1) to compute the multivariate 
rational generating function of C in polynomial time for varying dimension. 

To obtain a shorter formula, we use the technique of half-open exact decompositions [KV08], 
which is a refinement of the method of "irrational" perturbations [BS07b, Kop07b] . We use 
the following result; see also Figure 2.5 and Figure 2.6. 

Lemma 39. 

(a) Let 



be a linear identity (with rational coefficients £i) of indicator functions of cones Ci C M", 
where the cones Ci are full- dimensional for i £ Ii and lower- dimensional for i £ I2. Let 
each cone be given as 



(2.5) 




ie/2 



(2.6) 



= { X G : (b*^.,x) < forj G J, }. 



Let y G be a vector such that {h*j,y 
define the "half-open cone" 



) 7^ for all i £ IiU L2, j G J?,. For i G Ii, we 



(2.7) 




) <0, 




Then 



(2.8) 



ie/i 
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(modulo lower-dimensional cones) 



Figure 2.5. An identity, valid modulo lower-dimensional cones, correspond- 
ing to a polyhedral subdivision of a cone 




Figure 2.6. The technique of half-open exact decomposition. The relative 
location of the vector y (represented by a dot) determines which defining 
inequalities are strict (broken lines) and which are weak (solid lines). 



(b) In particular, let 

(2.9) [C] = + ^ [C] 

be the identity corresponding to a triangulation of the cone C, where T is the set of 
simplicial cones of maximal dimension and T is the set of intersecting proper faces. 
Then there exists a polynomial-time algorithm to construct a vector y € Q" such that 
the above construction yields the identity 

(2.10) [C]=5][C.], 

c,er 

which describes a partition of C into half-open cones of maximal dimension. 



Proof. Part (a) is a slightly less general form of Theorem 3 in [KV08]. Part (b) follows 
from the discussion in section 2 of [KV08] . □ 
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Since the cones in a triangulation T of all tangent cones C-p(v) of our polytopes are 
unimodular by Lemma 26, we can efficiently write the multivariate generating functions of 
their half-open counterparts. 

Lemma 40 (Lemma 9 in [KV08]). Let C C be an N- dimensional half- open pointed 
simplicial affine cone with an integral apex v G and the ray description 

(2.11) C = I V + Xlj^i Ajbj : \j>0 for j € J< and \j > for j S J< } 

where J< U J< = {1, . . . , N} and hj G Z" \ {0}. We further assume that C is unimodular, 
i.e., the vectors hj form a basis of the lattice (Mbi + • • • + Mbjv) H Z". Then the unique 
point in the fundamental parallelepiped of the half- open cone C is 

(2.12) a = v+J^b,-, 

ieJ< 

and the generating function of C is given by 

(2.13) gciz] 



Taking all results together, we obtain: 

Corollary 41. Let r be a fixed integer. There exist algorithms that, given 

(a) a matroid M on n elements, presented by an evaluation oracle for its rank function (p, 
which is bounded above by r, or 

(b) an evaluation oracle for an integral polymatroid rank function ip: 2["] — > N, which is 
bounded above by r, 

compute in time polynomial in n vectors aj G Z", bjj G Z" \ {0}, and Vj G for i £ I 
(a polynomial-size index set) and j = 1,...,N, where N < n, such that the multivari- 
ate generating function ofV{M), V^{M) and V{ip), respectively, is the sum of rational 
functions 

(2-14) gp(z) = y ;v ^ ' u — 
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and the k-th dilation of the polytope has the multivariate rational generating function 

2,ai + (fc-l)vi 



(2.15) gkvi^) = 



^nf=i(i--^-r 

Proof. Lemma 23 implies that finding the multivariate generating function of V{M), 
can be reduced to finding the multivariate generating functions of their 
tangent cones. Moreover, V{M), 'P-^(M) and V{ip) have only polynomially in n many 
vertices as described in subsection 2.3 or Lemma 35. Enumerating their vertices can be 
done in polynomial time by Lemma 17 or 36. 

Given a vertex v of P(M), P^(M) or P(V'), its neighbors can be computed in polynomial 
time by Lemma 17 or 36. The tangent cone C(v) at v is generated by elements in R^, where 
Ra is an elementary set for some A CI [n]. See subsection 2.3 or Lemma 37. We also proved 
in Lemma 26 that every triangulation of C(v) generated by elements in R^ is unimodular 
and Lemma 28 states that any triangulation of C(v) has at most a polynomial in n number 
of top-dimensional simplices. Moreover, a triangulation of the cone C(v) can be computed 
in polynomial time by Lemma 31. Finally, using Lemmas 39 and 40 we can write the 
polynomial sized multivariate generating function of C(v) in polynomial time. Therefore 
we can write the multivariate generating function of V{M), 'P-^(M) or V{ip) in polynomial 
time. □ 



2.6. Polynomial-time specialization of rational generating functions in varying 

dimension 

We now compute the Ehrhart polynomial i(V, k) = #{kV H Z") from the multivariate 
rational generating function g^-pi^) of Corollary 41. This amounts to the problem of eval- 
uating or specializing a rational generating function g^p^z), depending on a parameter k, 
at the point z = 1. This is a pole of each of its summands but a regular point (removable 
singularity) of the function itself. From now on we call this the specialization problem. We 
explain a very general procedure to solve it which we hope will allow future applications. 

To this end, let the generating function of a polytope V C M" be given in the form 

(2-16) 9v{z) = Y,e.,—-^^-— 

f^i njLi(l-z^'0 
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where et G {il}, G and hij E \ {0}. Let s = maxjg/ Si be the maximum number 
of binomials in the denominators. In general, if s is allowed to grow, more poles need to be 
considered for each summand, so the evaluation will need more computational effort. 

In previous literature, the specialization problem has been considered, but not in suffi- 
cient generality for our purpose. In the original paper by (author?) [Bar94, Lemma 4.3], the 
dimension n is fixed, and each summand has exactly Si = n binomials in the denominator. 
The same restriction can be found in the survey by (author?) [BP99] . In the more general 
algorithmic theory of monomial substitutions developed by (author?) [BW03, Woo04], 
there is no assumption on the dimension n, but the number s of binomials in the denomi- 
nators is fixed. The same restriction appears in the paper by (author?) [VW08, Lemma 
2.15]. In a recent paper, (author?) [Bar06, section 5] gives a polynomial-time algorithm 
for the specialization problem for rational functions of the form 



where the dimension n is fixed, the number s of different binomials in each denominator 
equals n, but the multiplicity 7ij is varying. 

We will show that the technique from (author?) [Bar06, section 5] can be implemented 
in a way such that we obtain a polynomial-time algorithm even for the case of a general 
formula (2.16), when the dimension and the number of binomials are allowed to grow. 

Theorem 42 (Polynomial-time specialization [DLIIK08]). (a) There exists an algorithm 
for computing the specialization of a rational function of the form 



at its removable singularity z = 1, which runs in time polynomial in the encoding size 
of its data Si G Q, a^ G for i & I and hij £ Z" for i £ I, j = 1, . . . , Sj, even when 
the dimension n and the numbers Si of terms in the denominators are not fixed. 
(b) In particular, there exists a polynomial-time algorithm that, given data £{ G Q, a-j G Z'^ 
for i £ I and hij G for i G /, j = 1, ... ,5, describing a rational function in the 
form (2.18), computes a vector X G Q" with (A, hij) / for all i,j and rational weights 




■ 3-1 



(2.17) 



(2.18) 
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Wi^l for i G I and I = 0, . . . , Si. Then the number of integer points is given by 

Si 

(2.19) #{VnZ^) = ^ei^Wi,i{X,Siiy. 

i£l 1=0 

(c) Likewise, given a parametric rational function for the dilations of an integral polytope V , 

2,ai+(fe-l)vi 



(2.20) g,^(^^) = Y^e.^^^ 

iai iij=i\^ 

the Ehrhart polynomial i(V, k) = ij^{kV H Z") is given by the explicit formula 

(2.21) i{V, fe) = E f E ^*)" E - 

m=0 ViG/ l=m ^"^^ / 

where M = min{s, dim'P}. 

Proof of Theorem 22. Corollary 41 and Theorem 42 imply Theorem (Theorem 22) 
directly. □ 

The remainder of this section contains the proof of Theorem 42. We follow [BP99] 
and recall the definition of Todd polynomials. We will prove that they can be efficiently 
evaluated in rational arithmetic. 

Definition 43. We consider the function 

a function that is analytic in a neighborhood of 0. The m-th (s-variate) Todd polynomial 
is the coefficient of in the Taylor expansion 

oo 

H{x,iu--.,is) = Y.^^rn{iu.-.,is)x'^. 

m=0 

We remark that, when the numbers s and m are allowed to vary, the Todd polynomials 
have an exponential number of monomials. 

Theorem 44 ([DLHK08]). The Todd polynomial tdm{S,i, ■■■ ,(s) can be evaluated for 
given rational data ^i, . . . , in time polynomial in s, m, and the encoding length o/^i, ■ ■ ■ ,^s- 

The proof makes use of the following lemma. 
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Lemma 45 ([DLIIK08]). The function h{x) = x/{\ — e~^) is a function that is analytic 
in a neighborhood ofO. Its Taylor series about x = is of the form 

oo 

(2.22) h{x) = V bnx'^ where bn = — —rrSn 

n!(n-|-l)! 

n=0 ' 

with integer coefficients Cn that have a binary encoding length o/0(n^ log n) . The coefficients 
Cn can be computed from the recursion 

Co = 1 

^'■'^^ Cn = X:(-l)^+^f^^) V "'\,,, cn-, for n = 1,2,.... 

Proof. The reciprocal function h^^{x) = (1 — e~^)/x has the Taylor series 

/i~-^(x) = V'onX" with On = 7 —- . 

^-^ (n + 1 

Using the identity h^^{x)h{x) = (X]^o '^n^") (Z^^o ^n^") = 1' obtain the recursion 

6o = - = 1 

(2.24) 

&n = -(oifen-i + a2^n-2 H V CLnbo) for n = 1, 2, 

We prove (2.22) inductively. Clearly bo = cq = 1. For n = 1, 2, . . . , we have 

Cn = n! (n + 1)! 6„ 

= -n! (n + 1)! (ai6„_i + a26„_2 H h a„6o) 

= + 1)! E J-TT)\ ■ (n-j)!(n-j + l)!'"-^' 

frC U + 1)1 {n-jy. {n-j + iy. 

Thus we obtain the recursion formula (2.23), which also shows that all Cn are integers. A 
rough estimate shows that 

\cn\ < n{n + ly. n\ |c„_i| < ((n + 1)!)^ |c„_i| , 

thus \cn\ < ((n + 1)!) ' , so Cn has a binary encoding length of 0(72^ log??.) . CH 
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Proof of Theorem 44. By definition, we have 

oo s 
m=0 j=l 

From Lemma 45 we have 

(2.25) h{xCj) = (^J,nx" + o(x'") where = -j^^—rr^Cn 

n=o n.[n + L). 

with integers c„ given by the recursion (2.23). Thus we can evaluate tdm(Ci, • • • , ■^s) by 
summing over all the possible compositions ni + • • • + = m of the order m from the 
orders Uj of the factors: 

(2.26) td„(6,...,es) = Yl Pl.n, . ■ ■ Ps,ns 

(ni,...,n^)eZ^ 
niH \'ns=m 

We remark that the length of the above sum is equal to the number of compositions of m 
into s non- negative parts, 

{m + s — l)(m + s — 2) . . . (m + s — (s — 1)) 
~ (s-l)(s-2)...2-l 

=^((1+5^; 



which is exponential in s (whenever m > s). Thus we cannot evaluate the formula (2.26) 
efficiently when s is allowed to grow. 

However, we show that we can evaluate tdm('Ci) •••)£«) more efficiently. To this end, 
we multiply up the s truncated Taylor series (2.25), one factor at a time, truncating after 
order m. Let us denote 

Hi{x) = h{xii) 
H2{x) = Hi{x) ■ h{xi2) 



Hs{x) = Hs-i{x) • h{x^s) = H{x,^i, . . . ,£s)- 
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Each multiplication can be implemented in O(m^) elementary rational operations. We 
finally show that all numbers appearing in the calculations have polynomial encoding size. 
Let H be the largest binary encoding size of any of the rational numbers , . . . , • Then 
every j3j^n given by (2.25) has a binary encoding size 0(Hn^ log'^ n). Let Hj{x) have the 
truncated Taylor series Y2^=o ^j,nx'^ + o(x'") and let Aj denote the largest binary encoding 
length of any aj^n for n < m. Then 

m n 

Hj+i{x) = '^aj+i^nx"' + o{x"^) with aj+i,n = y^^aj,i(3j^n-i- 

n=0 Z=0 

Thus the binary encoding size of aj+i,n (for n < m) is bounded by Aj + 0(Hm^ log^ m). 
Thus, after s multiplication steps, the encoding size of the coefficients is bounded by 
0(,sHm^ log'^ m), a polynomial quantity. □ 

Proof of Theorem 42. Parts (a) and (b). We recall the technique of (author?) 
[Bar94, Lemma 4.3], refined by (author?) [Bar06, section 5]. 

We first construct a rational vector A G such that (A,bjj) ^ for all One such 
construction is to consider the moment curve A(^) = (1, ^, . . . , ^"~^) G M". For each 
exponent vector bjj occuring in a denominator of (2.16), the function fij: ^ (A(^),bjj) 
is a polynomial function of degree at most n — 1. Since bjj ^ 0, the function fij is not 
identically zero. Hence fij has at most n — 1 zeros. By evaluating all functions fij for i £ I 
and j = 1, Si at M = (n — l)s|/| + 1 different values for for instance at the integers 
,^ = 0, . . . , M, we can find one ^ = ^ that is not a zero of any fij. Clearly this search can be 
implemented in polynomial time, even when the dimension n and the number s of terms in 
the denominators are not fixed. We set A = A(^). 

For T > 0, let us consider the points = e"^"^ = (exp{rAi}, . . . ,exp{rAn}). We have 



Zr'' = JJexp{rAi6ijz} = exp{r {X,hij)}; 
1=1 



since (A, b^j) ^ for all all the denominators 1 — z^'^ are nonzero. Hence for every 
r > 0, the point Zr is a regular point not only of g{z) but also of the individual summands 
of (2.16). We have 



g{l) = lim Vei 



^-'^^i n-Li(i-z>) 
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= lim exp{T(A,ai)} 
^-0+^ 'n;Li(l-exp{r(A,b.,)}) 

T 

= lim > r~*' expjr (A,aj)| TT — —— — 

--0+^ ^ ' *^a^ii_exp{r(A,b,,)} 

-1 

= lim Veir-"»exp{r(A,ai)}n — -— /i(-r(A,bij-)) 
^-0+ ^^^^^ (A, bi,-) 

= lim T-r«. /X u \ ^"'^ «^XP{^ a,)}/7(r, - (A, b^i) , . . . , - (A, b,,J) 

■^^0+-^ llj=l DijV 

where H(x,S^i, . . . is the function from Definition 43. We will compute the limit by 
finding the constant term of the Laurent expansion of each summand about r = 0. Now 
the function r ^ exp{r (A,aj)} is holomorphic and has the Taylor series 

(2.27) exp{T(A,a,)} = J^ai,y + o(t^0 where a,,; = 
and H{t, ■ ■ ■ ,S,Si) has the Taylor series 

/^(T,ei,...,6) = x;td'„(ei,...,6)T"^ + o(T^o- 

m=0 

Because of the factor r~*', which gives rise to a pole of order Sj in the summand, we can 
compute the constant term of the Laurent expansion by summing over all the possible 
compositions Si = I + {si — I) oi the order sf. 

(2.28) g{i) = E^^fpr^rvrr E ^^i) b,,J). 

We use the notation 

^ .Ns. td,,„/(-(A,b,,i),...,-(A,b,,,J) r ■ r .i_r. 

^'•'-^ /!.(A,b,,)...(A,b,,) for.e/and/-0,...,.., 

these rational numbers can be computed in polynomial time using Theorem 44. We now 
obtain the formula of the claim, 

Si 

5(1) = E E • 

iei 1=0 
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Part (c). Applying the same technique to the parametric rational function (2.20), we obtain 

Si 

= ^ ei ^ Wi,i (A, a.i + {k- l)vj)' , 



is/ 1=0 

Si I 



ig/ i=0 m=0 ^ ^ 

^ ( (A, V,)- f; ( M w,, (A, a. - v.)'-") 



m=0 l=m 

an explicit formula for the Ehrhart polynomial. We remark that, since the Ehrhart poly- 
nomial is of degree equal to the dimension of "P, all coefficients of /c™ for m > dimP must 
vanish. Thus we obtain the formula of the claim, where we sum only up to min{s,dim'P} 
instead of s. □ 
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CHAPTER 3 

Results on the Algebraic Combinatorics of Matroid 

Polytopes 

3.1. Algebraic Properties of /i*-vectors and Ehrhart polynomials of Matroid 

Polytopes 

Here we investigate algebraic properties of the Ehrhart functions of matroid polytopes: 
The Ehrhart series of a polytope V is the infinite series X^fc^o *('^' ^)^^- recall the 
following classic result about Ehrhart series (see e.g., [Hib92, Sta96]). Let V C M" be an 
integral convex polytope of dimension d. Then it is known that its Ehrhart series is a 
rational function of the form 



The numerator is often called the h* -polynomial of V (some authors also call it the Ehrhart 
h-polynomial) , and we define the coefficients of the polynomial in the numerator of Lemma 



3.1, hl + hlt-\ h /irf^i*'^"-^ + h^f^, as the h*-vector of V, which we write as h*{V) := 



A vector (cq, . . . , Cd) is unimodal if there exists an index p, < p < d, such that Cj_i < Cj 
for i < p and Cj > Cj^i for j > p. Due to its algebraic implications, several authors have 
studied the unimodality of /i*-vectors (see [Hib92] and [Sta96] and references therein). 

Suppose, as before, that V C M", and each vertex of 7^ has integral (or rational) vertices. 
Let Yi,Y2, . . . ,Yn and T be indeterminates over a field K. Letting g > 1 we define A(V)q 
as the vector space over K which is spanned by the monomials Y""^,!^"^, . . . ,Y^"T'^ such 
that (ai,a2; • • • ,««) ^ qP r\ Z". Since V is convex it follows that A{V)qA{V)p C Aq^p{V) 
for all p,q, and thus the Ehrhart ring of V, A{V) = ^ graded algebra 



(3.1) 




{hl,h\,...,h*,^,,hl). 



[Hib92, Sta96]. 
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It is well-known that if the Ehrhart ring of an integral polytope V, A{V)^ is Gorenstein, 
then h*i(P) is unimodal, and symmetric [Hib92, Sta96]. Nevertheless, the vector h*(P) can 
be unimodal even when the Ehrhart ring Ai(P) is not Gorenstein. For matroid polytopes, 
their Ehrhart ring is indeed often not Gorenstein. For instance, De Negri and Hibi [DNH97] 
prove explicitly when the Ehrhart ring of a uniform matroid polytope is Gorenstein or not. 
Two fascinating facts, uncovered through experimentation, are that all /i*-vectors seen 
thus far are unimodal, even for the cases when their Ehrhart rings are not Gorenstein. 
In addition, when we computed the explicit Ehrhart polynomials of matroid polytopes we 
observe their coefficients are always positive. We conjecture: 

Conjecture 46. Let V{M) he the matroid polytope of a matroid M. 

(A) The h* -vector ofV{M) is unimodal. 

(B) The coefficients of the Ehrhart polynomial ofV{M) are positive. 

We have proved both parts of this conjecture in many instances. A class of matroids 
that we considered are the uniform matroids; recall that the uniform matroid on n elements 
of rank r is the collection of all r-subsets of n. Using computers, we were able to verify 
Conjecture 46 for all uniform matroids up to 75 elements as well as for a wide variety of 
non-uniform matroids which are collected at [HawOD] . We include here this information just 
for the 28 famous matroids presented in [0x192]. Results in [Kat05], with some additional 
careful calculations, imply that Conjecture 46 is true for all rank 2 uniform matroids. 
Regarding part (A) of the conjecture we were also able to prove partial unimodality for 
uniform matroids of rank 3. Concretely we obtained: 

Theorem 47. 

(1) Conjecture 46 is true for all uniform matroids up to 75 elements and all uniform 
matroids of rank 2. It is also true for all matroids listed in [Haw09]. 

(2) Let ViU^'"") he the matroid polytope of a uniform matroid of rank 3 on n elements, 
and let I he a non-negative integer. Then there exists n{I) € N such that for all 
n > n{I) the h* -vector of V{U^'"'), (/iq, • • • , h'^), is non- decreasing from index to 
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Using the programs cdd+ [Fuk06], LattE [DLHH+05] and LattE macchiato [Kop07a] 
we explored patterns for the Ehrhart polynomials of matroid polytopes. Since previous au- 
thors proposed other invariants of a matroid (e.g., the Tutte polynomials and the invariants 
of [Spe06, BJR06]) we wished to know how well does the Ehrhart polynomial distinguish 
non-isomorphic matroids. It is natural to compare it with other known invariants. Some 
straightforward properties are immediately evident. For example, the Ehrhart polynomial 
of a matroid and that of its dual are equal. Also the Ehrhart polynomial of a direct sum of 
matroids is the product of their Ehrhart polynomials. 

Table 3.1. Coefficients of the Ehrhart polynomials of selected matroids in 
[BJR06, Spe06] 



Ehrhart Polynomial 



Speyerl 


1, 


21 

5 ' 


343 63 91 77 29 
45 ' 8 ' 18 ' 40 ' 90 


Speyer2 


1, 


135 
28 


3691 1511 88 39 529 89 
' 360 ' 120 ' 9 ' 8 ' 360 ' 420 


BJRl 


1, 


109 
30 


23 59 9 9 
' 4 ' 12' 4' 20 


BJR2 


1, 


211 
60 


125 33 43 43 
' 24 ' 8 ' 24' 120 


BJR3 


1, 


83 


2783 199 391 247 61 


20' 


360 ' 24 ' 72 ' 120' 180 


BJR4 


1, 


25 
6 ' 


353 101 193 23 53 
45 ' 12 ' 36 ' 12' 180 



We call the last two matroids in Figure 2 in [Spe06] Speyerl and Speyer2 and the 
matroids of Figure 2 in [BJR06] BJRl, BJR2, BJR3, and BJR4, and list their Ehrhart 
polynomials in Table 3.1. We note that BJR3 and BJR4 have the same Tutte polynomial, 
yet their Ehrhart polynomials are different. This proves that the Ehrhart polynomial cannot 
be computed using deletion and contractions, as is the case for the Tutte polynomial. 
Examples BJRl and BJR2 show that the Ehrhart polynomial of a matroid may help to 
distinguish non-isomorphic matroids: These two matroids are not isomorphic yet they have 
the same Tutte polynomials and the same quasi-symmetric function studied in [BJR06]. 
Although they share some properties, there does not seem to be an obvious relation to 
Speyer's univariate polynomials introduced in [Spe06]; examples Speyerl and Speyer2 show 
they are relatively prime with their corresponding Ehrhart polynomials. 

Our experiments included, among others, many examples coming from small graphical 
matroids, random realizable matroids over fields of small positive characteristic, and the 
classical examples listed in the Appendix of [0x192] for which we list the results in Table 
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3.2. For a comprehensive list of all our calculations visit [Haw09]. Soon it became evident 

that all instances verified both parts of our Conjecture 46. 

Table 3.2. Coefficients of the Ehrhart polynomials and h* -vectors of se- 
lected matroids in [0x192] 





n 


T 


/i*-vector 


Ehrhart Polynomial 




6 


3 


1, 10, 20, 10, 1 


1, 


107 21 49 7 7 
30 ' 4 ' 12 ' 4 ' 20 




6 


3 


1, 11, 24, 11, 10 


1, 


18 11 992 
5 ' 2 ' 2 ' ' 5 


Qe 


6 


3 


1, 12, 28, 12, 1 


1, 


109 23 59 9 9 
30 ' 4 ' 12 ' 4 ' 20 


Pe 


6 


3 


1, 13, 32, 13, 1 


1, 


11 R 16 5 1 
3 ' ' 3 ' 2 ' 2 


Re 


6 


3 


1, 12, 28, 12, 1 


1, 


109 23 59 9 9 
30 ' 4 ' 12 ' 4 ' 20 


F7 


7 


3 


21, 98, 91, 21, 1 


1, 


21 343 63 91 77 29 
5 ' 45 ' 8 ' 18 ' 40 ' 90 


P7 


7 


3 


21, 101, 97, 22, 1 


1, 


253 2809 33 193 61 121 
60 ' 360 ' 4 ' 36 ' 30 ' 360 




7 


3 


21, 104, 103, 23, 1 


1, 


127 479 69 17 257 7 
30 ' 60 ' 8 ' 3 ' 120 ' 20 


AG(3, 2) 


8 


4 


1, 62, 561, 1014,449, 48, 1 


1, 


209 1981 881 119 95 499 89 
42 ' 180 ' 60 ' 9 ' 12 ' 180 ' 210 


AG' (3,2) 


8 


4 


1, 62, 562, 1023, 458, 49, 1 


1, 


299 4007 5401 122 2911 1013 77 
60 ' 360 ' 360 ' 9 ' 360 ' 360 ' 180 




8 


4 


1,62,563,1032,467,50,1 


1, 


524 1013 1379 125 743 257 136 
105 ' 90 ' 90 ' 9 ' 90 ' 90 ' 315 


Fs 


8 


4 


1,62,563,1032,467,50,1 


1, 


524 1013 1379 125 743 257 136 
105 ' 90 ' 90 ' 9 ' 90 ' 90 ' 315 


Qs 


8 


4 


1, 62, 564, 1041, 476, 51, 1 


1, 


2099 4097 1877 128 337 1043 61 
420 ' 360 ' 120 ' 9 ' 40 ' 360 ' 140 


Ss 


8 


4 


1, 44, 337, 612, 305,40, 1 


1, 


1021 377 475 193 511 65 67 
210 ' 36 ' 36 ' 18 ' 90 ' 36 ' 252 


Vs 


8 


4 


1,62,570,1095,530,57,1 


1, 


2117 4367 2107 146 1133 1133 193 
420 ' 360 ' 120 ' 9 ' 120 ' 360 ' 420 


Tg 


8 


4 


1, 62, 564, 1041, 476, 51, 1 


1, 


2099 4097 1877 128 337 1043 61 
420 ' 360 ' 120 ' 9 ' 40 ' 360 ' 140 




8 


4 


1,62,569,1086,521,56,1 


1, 


151 2161 3103 143 1669 559 41 
30 ' 180 ' 180 ' 9 ' 180 ' 180 ' 90 


Ls 


8 


4 


1,62,567,1068,503,54,1 


1, 


527 529 83 137 134 136 47 
105 ' 45 ' 5 ' 9 ' 15 ' 45 ' 105 


J 


8 


4 


1,44, 339, 630, 323,42, 1 


1, 


512 193 83 205 361 17 23 
105 ' 18 ' 6 ' 18 ' 60 ' 9 ' 84 


Pa 


8 


4 


1, 62, 565, 1050, 485, 52, 1 


1, 


1051 2071 2873 131 1547 529 277 
210 ' 180 ' 180 ' 9 ' 180 ' 180 ' 630 


Wi 


8 


4 


1, 38, 262,475, 254, 37, 1 


1, 


135 3691 1511 88 39 529 89 
28 ' 360 ' 120 ' 9 ' 8 ' 360 ' 420 




8 


4 


1, 38, 263,484, 263, 38, 1 


1, 


169 467 581 91 227 68 68 
35 ' 45 ' 45 ' 9 ' 45 ' 45 ' 315 




9 


5 


78, 1116, 3492, 3237, 927, 72, 1 


1, 


307 137141 3223 37807 211 5743 1889 8923 
56 ' 10080 ' 160 ' 1920 ' 16 ' 960 ' 1120' 40320 


AG(2,3) 


9 


3 


1, 147, 1230, 1885, 714, 63, 1 


1, 


1453 41749 581 34069 927 4541 239 449 
280 ' 3360 ' 32 ' 1920 ' 80 ' 960 ' 224' 4480 


Pappus 


9 


3 


1, 147, 1230, 1915, 744, 66, 1 


1, 


729 3573 381 1499 243 49 153 57 
140 ' 280 ' 20 ' 80 ' 20 ' 10 ' 140 ' 560 


Non-Pappus 


9 


3 


1, 147, 1230, 1925, 754, 67, 1 


1, 


4379 25951 9287 21967 987 2855 3701 275 
840 ' 2016 ' 480 ' 1152 ' SO ' 576 ' 3360' 2688 


Q3(GF(3)*) 


9 


3 


1,147, 1098,1638,632,59,1 


1, 


433 3079 4193 5947 167 601 787 149 
84 ' 252 ' 240 ' 360 ' 16 ' 144 ' 840 ' 1680 


R9 


9 


3 


1, 147, 1142, 1717, 656, 60, 1 


1, 


723 49 88 24217 1291 625 821 133 
140 ' 4 ' 5 ' 1440 ' 120 ' 144 ' 840 ' 1440 



By far the most comprehensive study we made was for the family of uniform matroids. 
In this case we based our computations on the theory of Veronese algebras as developed 
by M. Katzman in [Kat05]. There, Katzman gives an explicit equation for the /i*-vector of 
uniform matroid polytopes (again, using the language of Veronese algebras). For this family 
we were able to verify computationally the conjecture is true for all uniform matroids up 
to 75 elements and to prove partial unimodality as explained in the introduction. 

Lemma 48. The coefficients of the Ehrhart polynomial of the matroid polytope of the 
uniform matroid [7^'" are positive. 
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Proof. We begin with the expression in Corollary 2.2 in [Kat05] which explicitly gives 
the Ehrhart polynomial of ({/'''") as 

(3.2) m,r-u) =£(-!)• (;;) ('•<'' - • + 

Letting r = 2, Equation (3.2) becomes 

(2k + n- l)(2k + n - 2) ■ ■ ■ (2k + 1) (k + n - 2)(k + n - 3) ■ ■ ■ (k) 
(n-l)! (n-1)! 

We next consider the coefficient of A;"""^"^ for < p < n — 1, which can be written as 

2„-p-i y- 1 

^-^ (n — 7i ) • • • (ra — iri) 



n 



E 



n — 1 ^ (n — 7i) • • • (n — i-n-x) 

y ( ^ Y ^ 

Kin — ji) ■ ■ ■ (n — j-o-i) ^ n — j„ 

n V - 1 

n — 1 (n — 7i) ■ ■ ■ (n — 7b-i) 



(3-3) = E ( 7 2"-^'-! V — 



n-l 

n 



n — jr, n — 1 

jp=i+jp-i 



It is known that the constant in any Ehrhart polynomial is 1 [Sta96], thus we only need 
to show that Equation (3.3) is positive for < p < n — 2. It is sufficient to show that the 
square- bracketed term of (3.3), 

n-l 



(3.4) 2"-P-i V — 



1+ ' 



n — jr, V n — 1 

is positive for < p < n — 2 and all 1 < ji < ■ ■ ■ < jp^i < n — 2. We can see that 
y^"7L\ I — ^ > 1. Moreover, since < p < n — 2 then 2"~p~^ > 2 and 1 H — < 2 since 

^-'Jp— -"-"TJp— 1 n— — ' — — — n—L — 

n > 2, proving the result. □ 
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To present our results about the /i*-vector we begin explaining the details with the 
following numbers introduced in [Kat05], which we refer to as the Katzman coefficients: 

Definition 49. For any positive integers n and r define the coefficients A^'^ by 

n(r— 1) 

^^''^r = (i + r + --- + r'-i)". 

We also define the vector A"'*" as (^A^'^ , ^4"''", . . . , A^'^^^-^^^ . 

Looking at the definition of the Katzman coefficient, we see that Aj'^ = and A"'^ = 
unless j = 0, in which case we have A^'^ = 1. 

Below we derive some new and useful equalities and prove symmetry and unimodality 
for the Katzman coefficients. Katzman [Kat05] gave an explicit equation for the /i*-vector 
of uniform matroid polytopes and the coefficients of their Ehrhart polynomials, although 
he did not use the same language. We restate it here for our purposes: 

Lemma 50 (See Corollary 2.9 in [Kat05]). Let P(C/^'") be the matroid polytope of the 
uniform matroid of rank r on n elements. Then the h* -polynomial of V{U^''^) is given by 



r-l 



(3.5) 



s=0 j=0 A:=0 l>k 



That is, for h*(P(C/"''^))) = {h*^, ...,h*), 

r-l s j 



hi = j2Yl E(-i)'"''^' A""-'''- 



3=0 j=0 k=0 

For r = 2 the h* -polynomial of V{U)^''^ is 



s)[j)\kj'^(i-'^)(--^y 



(3.6) 



E 

l>0 



V.t' I - nr. 

21 



The following lemma is a direct consequence of Corollary 2.9 in [Kat05]: 



Lemma 51. Let ^([7^'") be the matroid polytope of the uniform matroid of rank 2 on n 
elements. The h* -vector of V{U^''^) is unimodal. 
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The rank 2 case is an interesting example already. Although the /i*-vector is unimodal, 
it is not always symmetric. Next we present some useful lemmas, the first a combinatorial 
description of the Katzman coefficients. 



Lemma 52. For i = 0, . . . , n{r — 1) we have 



(3.7) AY = 



Oao + laiH \-{r—l)ar-i=i 

ao+aiH \-ar-\=n 



n 

ao, oi, • ■ ■ 5 Q^r-l 



where oq, . . • , o-r-i i^un through non-negative integers. 



Proof. Using the multinomial formula [Sta97] we have 



n(r— 1) 

J2 ^"'''r =(i + r + --- + r'^-i) 

1=0 



ao+iiH \-ar-i=n 



Sr^ I \'^aorpairp2a2 rp{r-l)ar-i 

^ \ ao, ai, • • • ) Olt-i) 



Sr^ ( " \2.0ao+lai+-+(r-l)ar_i 



fiO+ilH \-ar-\=n 

By grouping the powers of T we get equation (3.7). □ 



Next we present a generalization of a property of the binomial coefficients. The following 
lemma relates the Katzman coefficients to Katzman coefficients with one less element. 



Lemma 53. 
(3.8) AY= 

where we define Ap~^'^ := when p < or p > [n — l){r — 1) . 



fe=i— r+1 
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Proof. 



n{r~l) 
1=0 



(1 + T + • • • + T'^-i)" = (1 + T + • • • + + T + • • • + T^-i^ 



^(n-l)(r-l) 

^ ^n-i,r^i I (1 + T + . . . + r 



j=0 



(n-l)(r-l) {"-l){r-l) (n-l)(r-l) 

_ ^ j^-l,rj,i _|_ ^ ^n-l.rjij+l _(_..._(_ ^ ^"-l.^j^i+r-l 

j=0 i=0 i=0 

(n-l)(r-l) (n-l)(r-l)+l (n-l)(r-l)+r-l 

j=0 i=l i=r—l 
(ra-l)(r-l)+r-l 

Thus we get equation (3.8). □ 



The following lemma relates the Katzman coefficients of rank r with those of rank r — 1 . 



Lemma 54. 



n y ^ (k{r-2) \ 



or in other words 



^7= E [ M 



k+l=i 
0<k<n 
0<l<k{r-2) 



Proof. From Definition 49 



n(r— 1) 

A^T =(1 + T + • • • + T^-^Y = (1 + [T + • • • + T"-^] )" 

i=0 

=E(^)[^+---+^^-^]'=EQ^ni+---+^^ 

A:=0 ^ ^ A;=0 ^ ^ 



-21 



n . s /fc{r-2) 

E(>M E ^r'T^ 

k=0 ^ ^ V z=o 



□ 
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Lemma 55. The Katzman coefficients are unimodal and symmetric in the index i. That 
is, the vector (^Aq'^ , A^'^ , . . . , is unimodal and symmetric. 

Proof. We first prove symmetry. Considering equation (3.7) we assume that 

Oao + lai + • • • + (r — l)ar-i = i and ao + ai + • • • + ar-i = n. 

These two assumptions imply that 

((r - 1) - 0)ao + ((r - 1) - l)ai + • • • + ((r - 1) - (r - l))a,_i 

=(r - l)ao + (r - l)ai H h (r - l)a„ - Oao - loi (r - l)ar_i 

=(r — l)n — i. 

Therefore 



Oao+laiH h(r'— l)ar_i=i 

ao+^iH (-ar-^i=n. 



n 

ao, ai, • • • , Or-i 



(r-l)ao + (r-2)aiH |-Oar_i=i 

lO+QiH har-l='^ 



n(r— ■!)—«■ 



To prove unimodality we proceed by induction on n, where r is fixed. First, A^'*" is 
unimodal. Assume for n — 1 that A""-'^''' = (^Aq~^'^ , A^~^'^ , ■ ■ ■ , A'^~^'-^^^^_-^-^^ is unimodal. 
Using equation (3.8) and the fact that A""^"'" is symmetric we get that A"''' is unimodal. 
To help see this, one can view equation (3.8) as a sliding window over r elements of the 
vector A""-'^''', that is, A^'^ is equal to the sum of the r elements in a window over the 
vector A"""'^''^'. As the window slides up the vector A""^''', the sum will increase. When 
the window is on the center of A"^-"^'*" symmetry and unimodality of A""-'^''' will imply 
unimodality of A"'**. □ 

Now we use the explicit equation for the /i*-vector of uniform matroid polytopes to 
prove partial unimodality of rank 3 uniform matroids. First we note that the coefficient of 
T\ h*i, in equation (3.5) is 

r— 1 s j 



s=0 j=0 k=0 
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Letting the rank r = 3, and using equation (3.5), we get the /i*-polynomial (which is 
grouped by values of s from (3.5)), 



l>0 



n 



An,2 An-1,2 , 



Now using that ^"'^ = (^) and A^'^ = 6o{i), where (5j(p) 



1 if p = j 
else 



we get 



E 

l>0 



n\ j n — 1 
21 2/ 



n — 1 
21-2 



+ 



n 



+ ( 2 I - mi) + mi) + Soil) - 2S,{1) + 52(1)) 



E 

l>0 



n\ fn — 1 
2l) 21 



Using properties of the binomial coefficients, we see that 



+ 



n — 1 
21 



n — 1 
21-2 



n — 1 
21 

n — 1 
21 - 1 

n 

21 - 1 



n — 1 
21 - 1 

n — 1 
21-2 



+ 



n — 1 
21 



n — 1 
21-2 



So the /i* -polynomial of rank three uniform matroid polytopes is 



(3.9) 



E 

Z>0 



rpl 



Using Lemma 3.8, the coefficient of T\ if 3/ < n, is 



(3.10) 



hi 



E 



0<p<fc<n 



k J \j) 



21 - 1 



n| . 1 +(52(0 
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Next we show that when g is fixed A^'^ is a polynomial of degree g in the indeterminate 
n, with positive leading coefficient. Assume g < n. Considering Lemma 54 and when g <n, 



(3.11) 




0<p<k<n 



where (^) is a polynomial of degree q with positive leading coefficient. The highest degree 
polynomial in the sum is , a degree g polynomial. Hence A^''^ is a polynomial of degree 
g in the indeterminate n, with positive leading coefficient. If 5 > n, then AJg^ = A^'tg since 
the Katzman coefficients are symmetric by Lemma 55. 

Proof of Theorem 47 Part (2). Let / be a non-negative integer. From above we 
see that A^'^ is a degree g polynomial in the indeterminate n, with positive leading coef- 
ficient. Equation (3.9) is the /i* -polynomial of {7'^'", which is a sum of polynomials in n, 
the highest degree polynomial being ^4^^'^, a polynomial of degree 3/. So, /i^ — is the 
difference of a degree 3/ and 3(/ — 1) polynomial, hence — hj^__i is a degree 3/ polynomial 
with positive leading coefficient. For sufficiently large n, call it n(I), /i^ — h'^_i is positive for 
< / < /. Hence, the /i*-vector of C/^'" is non-decreasing up to the index / for n > n(/). □ 

One might ask if equation (3.11) has a simpler form. We ran the WZ algorithm on our 
expression, which proved that equation (3.11) can not be written as a linear combination of 
a fixed number of hypergeometric terms {closed form) [PWZ96]. There is still the possibility 
that this expression has a simpler form, though not a closed form as described above. 

3.2. Unimodular Simplices of Matroid Polytopes 

A k-simplex S = {si, . . . , s^+i} is the collection of -|- 1 affinely independent points in 
and we say conv(S') := conv(si, . . . ,Sfc+i). When we refer to a simplex S we consider 
it both as a collection of points and as a convex polytope when the context is clear. Let 
V C M" be a polytope with vertices vi, . . . , Vp. A triangulation T of is a collection of 
simplices on the vertices vi, . . . , Vp of P such that 



3.2. UNIMODULAR SIMPLICES OF MATROID POLYTOPES 



54 



(i) A £ T then all faces of A are in T, 

(ii) E A,B £T then An B is a face of both A and B, 

(iii) UA6TConv(A) =P. 

It is sufficient to list only the highest dimensional simplices of a triangulation. A trian- 
gulation T is unimodular if the volume of every highest dimensional simplex of T are the 
same. We also say a simplex is unimodular if its normalized volume is one [DLRS09]. A 
vector in i?" is a zizl-vector if it contains +1 and —1 and the rest are zeros. 

Recall that two elements s, t G S of a matroid Ai = {S, I) are connected if there exists a 
circuit that contains s and t. This defines an equivalence relation on S and moreover M can 
be written as a direct sum of restrictions to the connected components of M [0x192] . This 
implies a matroid polytope Vm can be written as a direct product of connected matroid 
polytopes. We note that if a polytope 7^ is a direct product of polytopes Qi, ■ ■ ■ , Qk which 
each have a unimodular triangulation, then V has a unimodular triangulation [DLRS09]. 
Hence, we need only consider connected matroids when investigating unimodular triangu- 
lations. 

We wish to find sufficient conditions such that a (n— l)-simplex of a connected matroid 
polytope Vm is unimodular. We will see later that a simplex of Vm is unimodular if and 
only if the determinant of the vertices is rank(A^). 

Neil White proposed an algebraic conjecture on the toric ideal determined by a matroid 
[WhiSO]. Let he a field and define the polynomial ring Sm '■= ^[yb \ B £ Bm]- Let the 
toric ideal Im he defined as the kernel of the homomorphism Sm ~^ J^[xi, • • ■ , Xn] given by 
Yb ^ HieB Xi. Given any pair of bases Bi, B2 £ B the symmetric exchange property states 
that for every x £ Bi there exists y £ B2 such that Bi — x + y and B2 — y + x are bases 
and it is said that x £ Bi double swaps into B2. If Bi — x + y and B2 — y + x are bases 
they are called a double swap. White proposed the following conjecture: 

Conjecture 56 ([WhiSO]). For any matroid M, the toric ideal Im is generated by the 
quadratic binomials yB^ysj ~ yD^yD2 such that the pair of bases Di,D2 can be obtained 
from the pair Bi, B2 by a double swap. 
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A stronger conjecture is that the previous quadratic binomials are in fact a Grobner 
bases for some term ordering. This imphes that there exists a regular unimodular triangu- 
lation of any matroid polytope. The following is a geometric variation of White's famous 
conjecture: 

Conjecture 57. Let M be a matroid of rank r on n elements. There exists a unimod- 
ular triangulation of Vm ■ 

In [Stu96] it was shown that all uniform matroids have a regular unimodular triangu- 
lation. A partial proof of White's original conjecture was given in [Bla08] for graphical 
matroids, but Conjecture 57 has yet to be proved in general. 

A covering of "P is a triangulation without condition (ii). That is, the simplices of a 
covering can intersect in any fashion, as long as the polytope is covered. See Figure 3.1. If 
asking for a unimodular triangulation is too much, we propose a weaker conjecture: 




Figure 3.1. A covering of a hexagon. 



Conjecture 58. Let M he a matroid of rank r on n elements. There exists a unimod- 
ular covering ofFju- 

Coverings are important to matroid polytopes. If a unimodular covering exists for a 
connected matroid polytope Vm it would imply its Caratheodory rank is n. In other words, 
every integer point in the cone of Vm can be written as an integral linear combination of 
vertices of Vm- That is, the vertices of Vm &re a Hilbert basis for the cone of Vm- See 
[CS86, Seb90, dPS03]. 
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We now focus on sufficient conditions for the existence unimodular simplices of matroid 
polytopes and begin by defining two graphs on any collection of incidence vectors of a 
matroid polytope. 

Definition 59. Let be a matroid on n elements. Given a subset X C Inc(;Bx) C 
we define two graphs G and g onn nodes. The graph G{X) has an edge {i, j) for every pair 
of bases Xj, Xj G X that are adjacent, i.e. Xj — Xj — — et for some s, t. The graph g{X) 
has an edge every pair of adjacent bases Xs,Xt G X such that they differ in position 
i and j, i.e. X^ — Xj = — e^. Succinctly we define G{X) and g{X) as: 

G{X) := ([n], { I Xi,X,- G X, X^ - X,- = e, - e* }) 

g{X) := {[n], { \ X^X* G X, X, - X* = - e,- }) 

where e^, 6^,65,64 G M"" are standard unit vectors. 

Throughout this section we will take X to mean a set of incidence vectors as well as 
a 0/1 matrix with rows as the incidence vectors. The following proposition is easy to see 
using elementary row sums on X and cofactor expansion on the resulting last row. 

Proposition 60. Let M be a connected matroid of rank r on n elements and X c 
Id.c{Bm) be n linearly independent incidence vectors oJVm- Then \det{X)\ = kr, A; G N. 

Example 61. Consider the following incidence vectors X in row form. We label the 
rows 1, . . . , 6 and columns a, . . . , f . 





a 


b 


c 


d 


e 


f 


1 


fl 


1 








1 




2 


1 


1 











1 


3 


1 





1 


1 








4 





1 


1 


1 








5 








1 





1 


1 


6 










1 


1 





G(X) = ([6],{(1,2),(3,4),(5,6)}) 
5(X) = ([6],{(a,b),(c,d),(e,f)}) 



We can equivalently define g{X) as the graph whose incidence matrix is given by 
{ Xg — Xt I Xg, Xj G X, X5 — Xj = — Gj } 
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where ej, ej G M" are standard unit vectors. We will use this fact for the following theorem. 

Theorem 62. Let M be a connected matroid of rank r onn elements and X c 1ivc{Bm) 
be n linearly independent incidence vectors of Bj^ . Then 

(i) G{X) and g{X) have the same number of connected components c, 

(ii) Let X = {v^, . . . , v"} and without loss of generality assume {1, . . . , h—l}, {Zi, . . . , I2— 

1), ■ ■ ■ , {Ic-i, ■ ■ ■ ^n}, are the connected components of g{X). Then 



(3.12) 



|det(X)| 



det 



El 



mi 



i = lr. 



.rn.2 



where v"^i , • . • , v"^= are representative incidence vectors of the connected compo- 
nents of G{X) . 



Proof. Wc offer a constructive proof of the theorem. Assume G{X) has p connected 
component and g{X) has c connected components. First let v"*i , . . . , v"**" be representa- 
tive incidence vectors of the connected components of G{X). Performing elementary row 
operations such as row addition, difference or reordering will only change the sign of the 
determinant. First we organize the rows of X by the p connected components of G{X). 



( < 



X 



Consider a connected component of G{X) with representative incidence vector v'"*. Let 
v"*' be the root of a spanning tree T of the corresponding component of G{X). For all non- 
root vectors in the component determined by v"** perform the row operation := — v'* 
where h is the parent of g. Now all the rows except v"** of the component of G{X) are 
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ibl-vectors. Repeat for each connected component of G{X). (By • . . i 1 • • • we mean some 
±1 vector). 



/m 
^1 



...±\... 
...±\... 



■ ±i- 



■ ±i- 




V 



■ ±i- 



■ ±i- 



Next, due to the assumption on the ordering of the connected components of giX) we 
partition the columns of matrix X by the components of g{X\ We also reorder the rows with 
the representatives vectors v™i, . . . , v™^ on top and the ±l-vectors by their corresponding 
components of g{X). 



mi mi 

"i •••"iiLi 

mo mo 

"i •••"ii-i 


mi m-i 

mo mn 

V, ■ ■ -v, ^ , 

'1 '2-1 

v"''' ■ ■ ■ t)!"" , 

'l '2-1 




"^1 mi \ 

V, • • • u' 2 

'c-1 

v7^^ ■ ■ ■ 

'c~l " 


...±i... 
...±i... 


0---0 
0---0 




0---0 
0---0 


0---0 
0---0 


...±1... 
...±1... 




0---0 
0---0 










0- • • 
0---0 


0---0 
0---0 




...±1... 



We note that for isolated nodes of g{X) there will be no corresponding ±l-vectors, but 
this will not interfere with this proof. Consider any submatrix of ibl-vectors determined 
by a connected component C of g{X). First, if there are A; many nodes in C, then there 
will be A; — 1 itl-vectors since X is linearly independent and these vectors are a spanning 
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tree for an orientation of the component C. Using the connectedness of the component C, 
elementary row addition and the fact that ±l-vectors sum to we can rewrite the parts of 
each representative incidence vector , • • • , contained in the columns corresponding 
to C. We sum the elements of each representative incidence vector over the component C 
and place it in the first position, with O's in all other position of C. 



■ ±i- 



..±1.. 

0---0 



0---0 



0---0 



0---0 



'2-1 
l2-l 



"'Til \ 



o-o 



0---0 

..±1.. 



■ ±1- 



0---0 



0---0 



0---0 



0---0 
0---0 



0---0 



■ ±1- 



± 1 



e!lt'C o-.-o 



±1- 



±1- 



0---0 



0---0 



V 



0---0 



0---0 



ElL^^^r o-.-o 



0...0 



0...0 



0. . . 



± 1. 



± 1. 



0. . . 



Er=,,_,''ro...o 

EIL,^,!".""' 0...0 



0...0 



0. . . 



0...0 



0...0 



. ± 1. 



■ ± 1. 



For each component C of g{X) there exists a column corresponding to a leaf of the spanning 
tree determined by the ±1 vectors in C, rooted at the first column. This column contains 
all O's and one 1. Performing cofactor expansion down this column eliminates a row and 
column. This can be repeated until all ztl-vectors of C are eliminated, preserving the 
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summed parts of the incidence vectors. Repeating for all connected components, all ±1- 
vectors can be eliminated while preserving the determinant of X (under a sign change). We 
get the desired form and since each operation preserves the matrix being square we see that 
c = p. 



\ 



□ 



Corollary 63. Let M he a connected matroid of rank r on n elements and X C 
Inc(i3_A4) he n linearly independent incidence vectors of Bjn- If G{X) is connected then 
\det{X)\ = r. 



Example 64. We offer an example of Theorem 62 in action. Consider the uniform 
matroid J7^'^ and the following collection of incidence vectors of bases in row form. 
Partition into components of G{X) and g{X). 



1 


1 








1 







1 


1 








1 





1 


1 











1 




1 


1 











1 


1 





1 


1 










1 





1 


1 











1 


1 


1 













1 


1 


1 














1 





1 


1 










1 





1 


1 











1 


1 


1 













1 


1 


1 



Find rooted spanning trees of each component of G{X) and take differences to get ±1- 
vectors and the representative vectors. 



1 


1 








1 







2 











1 

















-1 


1 
















-1 


1 


1 





1 


1 










1 





2 











-1 


1 
















-1 


1 




















1 





1 


1 










1 





2 











-1 


1 
















-1 


1 
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For each component of g{X) eliminate use itl-vectors to eliminate rows and columns using 
cofactor expansion. 



2 











1 

















-1 


1 


1 





2 











-1 


1 




















1 





2 











-1 


1 










2 








1 
















-1 


1 




1 


2 
















1 





2 










-1 


1 









2 





1 











-1 


1 


1 


2 











1 


2 






2 1 

0-11 

1 2 
10 2 
0-11 



2 10 
0-11 
12 
12 



2 1 
12 
12 



We explore more connections between the graphs G{X) and g{X). The following is a 
alternate proof of part of Theorem 62 

Lemma 65. Let M be a connected matroid of rank r on n elements and X C Inc(;B^) 
he n linearly independent incidence vectors of Bm- If G{X) is connected then g{X) is 
connected. 

Proof. Consider the matrix X and, as in the proof of Theorem 62, we pick a row 
as a root of the spanning tree of G{X) and perform row operations so that the resulting 
matrix is our root and the rest are itl-vectors. If g{X) is not connected then let S be one 
connected component and T = [n]\S. Also, let s = \S\ and t = \T\. Then the following 
non-zero vector is orthogonal to all the elements of X, contradicting its full- dimensionality. 



ie5 



□ 



More generally we can show the following: 



Lemma 66. Let M he a connected matroid of rank r on n elements and X C Inc(;S;v!) 
he n incidence vectors of M. If G{X) is connected, then rank(X) = n + 1 — ^connected 
components of g{X). 
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Proof. We will show that the space of functionals constant on X has dimension equal 
to the number k of connected components of g{X) minus one. Considering X as a matrix 
with incidence row vectors we first observe that all functionals in the kernel of X must give 
the same value to all entries (ej's) in the same connected component of g{X). That is, they 
are characterized by k parameters. Thus, the dimension of the space of functionals constant 
on X is at most k — 1. 

Consider the matrix X and, as in the proof of Theorem 62, we pick a row as a root 
of the spanning tree of G{X) and perform row operations so that the resulting matrix is 
our root root and the rest are ibl-vectors. Omitting the root row of X, the orthogonal 
complement of the remaining rows of X is exactly the A;-dimensional linear space described 
above. Reinserting the root row decreases the dimension of the orthogonal compliment by 
at most one. □ 

Lemma 67. Let M be a connected matroid of rank r on n elements and X C Inc(S^) 
be n incidence vectors of Bm- U di^) connected, rank(X) = n. Moreover G{X) is 
connected and thus \ det(X)| = r. 

Proof. Consider the matrix X and, as in the proof of Theorem 62, we pick representa- 
tive vectors v"^i , • • • , v'"*' for each of the components of G{X) . Next perform row operations 
so that the resulting matrix consists of the representative vectors and itl-vectors. Moreover 
since g{X) is connected these itl-vectors must be a spanning tree of an orientation of g{X), 
hence there must be (n — 1) itl-vectors. This implies that we have only one representative 
vector and thus G{X) is connected. We also see that these itl-vectors must be linearly 
independent. It follows that X must be linearly independent. Then |det(X)| = r follows 
from Corollary 63. □ 

In order to prove a unimodular covering exists we propose another conjecture which 
implies Conjecture 58: 

Conjecture 68. Let M be a matroid of rank r on n elements. For every x G Vm there 
exists n linearly independent vectors vi, . . . , v„ G Inc^Bjn) such that x G conv(vi, . . . , v„) 
and G{X) connected. Hence, conv(vi, . . . is unimodular. 
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In support of Conjecture 57, Conjecture 58 and Conjecture 68 we now present some 
experimental results. We obtained bases description for all matroids of nine elements or 
less from [MRar] . The number of matroids by rank and elements less than or equal to nine 
is given by the table: 



Elements 
Rank 





1 


2 


3 


4 


5 


6 


7 


8 


9 





1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 




1 


2 


3 


4 


5 


6 


7 


8 


9 


2 






1 


3 


7 


13 


23 


37 


58 


87 


Q 
O 








i 


4 


io 


oo 


iUo 




1 07^ 
IZ 1 


4 










1 


5 


23 


108 


940 


190214 


5 












1 


6 


37 


325 


190214 


6 














1 


7 


58 


1275 


7 
















1 


8 


87 


8 


















1 


9 


9 




















1 


Total Matroids 


1 


2 


4 


8 


17 


38 


98 


306 


1724 


383172 



With this data in hand we used TOPCOM [Ram02] to compute triangulations of all 
the connected matroid polytopes. From these computational experiments we can justify 
the following theorem: 

Theorem 69. All 1317 matroid polytopes of connected matroids with eight or less ele- 
ments have a unimodular triangulation. Moreover, the placing triangulation, determined by 
the order of bases as presented in [MRar], is a unimodular triangulation. 

Whether the ordering is important for the placing triangulation to be unimodular for 
these experiments is unclear. As a note, there certainly are non-unimodular triangulations 
of matroid polytopes. As a natural follow up experiment, we attempted to find as many 
matroids in the test set that were unimodular as well as triangulated by simplices that had 
connected graphs G{X). 
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Theorem 70. All 48 matroid polytopes of connected matroids with six or less elements 
have a triangulation where for each simplex X, G{X) is connected. Hence there exists a 
unimodular triangulation of each matroid polytope. 

For matroids with seven or more elements, it becomes especially time consuming to 
find a triangulation with G-connected simplices. We were able to find 36 matroids of rank 
seven that had G-connected simplices. All matroids and triangulations for Theorem 69 and 
Theorem 70 can be found at [Haw09]. 

3.3. Two Dimensional Faces of Matroid Polytopes 

In the pursuit of the proof of Conjecture 57, Conjecture 58, and Conjecture 68 related 
to matroid polytope triangulations, coverings and decompositions, we studied the two di- 
mensional faces of Vm order to understand how such subdivisions interact with the low 
dimensional faces. We obtained a few results. 

Lemma 71 (See Theorem 4.8 in [BGW07]). Let M he a matroid and Vm be its matroid 
polytope. The 2-dimensional faces ofVjn equilateral triangles or squares. 



Wi 




Figure 3.2. Square 2-face adjacency structure, e^ / e.^, e^ / e^, e^ ^ et, 

Lemma 72. Let M be a matroid of rank r on n elements and F be a 2-dimensional 
square face of Vm with vertices wi,W2,W3,W4 such that (wi,W2), (w2,W4), (wi,W3), 
and (w3,W4) are the edges. There exists standard unit vectors in es,et,em,ei G M" such 
that W2 = wi + e^ — e^, W3 = wi + e^ — e^, and W4 = wi + e^ — e^ + e^, — e/. Moreover, 
and Bs / e-m, et / e;, e^ / e^ and e^ / e^. 
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Proof. Let Ai he a matroid and F a 2-face of Vm with wi, W2, W3, W4 vertices of 
F such that (wi,W2), (w2,W4), (wi,W3), and (w3,W4) are the only adjacencies. Due to 
Lemma 17 we know if W2 = wi + — and W4 = wi + — + — e^, then W3 = 
wi + Bm — ei. Otherwise an additional adjacency (edge) is created among wi, W2, W3, W4, 
a contradiction to F a 2-face. E.g., if W3 = wi + — e/ then W3 and W2 are adjacent, a 
contradiction. Then by vector addition W4 = W3 +6^ — e/. Additionally, eg 7^ em, 7^ 
eg 7^ et and e,„ 7^ e^. Otherwise it would create an additional adjacency, hence edge. See 
Figure 3.2. 

□ 



/ • \ \ 
/ ' ' \ 

'W2 m ^ \ ")^W3 



Figure 3.3. Six adjacent vertices restrict any four from being a 2-face. 



Lemma 73. Let M be a matroid of rank r on n elements, wi,W2, W3,W4 G Vm- Letwi 
and W4 be adjacent to W2 and W3 in Vm where W2 = wi-|-(es— e^), W3 = wi-|-(em— e^) and 
W4 = wi (cs - et) + {em - ei), / e^, / e^, / et and em / e/. Then wi, W2, W3 
and W4 is a square 2-face if and only if'w^ := wi -|- (e^ — e^) or wg := wi -|- (e^ — e^) are 
not vertices ofVM- 

Proof. Assume that both W5 and wg are vertices of Vm and F is a square 2-face with 
vertices wi,W2,W3 and W4. Let 

(3.13) /(x) := aixi + 02X2 H h a„Xn = c 
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be a supporting hyperplane of F, i.e. wi,W2,W3,W4 satisfy (3.14). Moreover, for all 
z £ Vm \ f{z) < c. Since Equation (3.14) is a supporting hyperplane and W5, wg ^ F, 
/(ws) < c and /(we) < c. But 2c = /(w2 + W3) = /(w5 + we) < 2c, a contradiction. 



Assume, without loss of generality, that W5 is a vertex of Vmi while we is not. We will 
construct a supporting hyperplane of F. Let 

(3.14) /(x) := aixi + 02x2 H h a„x„. 

Let as = at = 1, am = a/ = 2, aj = 3 if i E supp(wi) \ {s,t,m,l} and aj = —1 if 
i i supp(wi) U {s, t, m, /}. We see that /(wi) = /(W2) = /(wg) = /(W4) = 3(r - 2) + 3, 
yet /(ws) = 3(r — 1) + 2 and for all vertices z G 'P_a4 \ /(z) < 3(r — 1) + 3. Hence / is a 
supporting hyperplane and F is a square 2-face. If W5 is also not a vertex, then / is still a 
supporting hyperplane of F . □ 

Corollary 74. Let M be a uniform matroid of rank r on n elements. Then the 
2-dimensional faces of Vm 0,^^ equilateral triangles. 

Proof. This follows from Lemma 73 and that if W2 = wi + (e^ — et) and W3 = 
wi + (em — e/) are vertices of Vm then so are W5 := v + (e^ — e^) and we := v + (e^ — e^) 
since is a uniform matroid. □ 

Proposition 75. Let M he a connected matroid of rank r on n elements and v he a 
vertex ofVM- Then all adjacent vertices of v lie on a {n — 2)-hyperplane 

Proof. Without loss of generality assume v = (1, . . . , 1, 0, . . . , 0)^. Then by Lemma 17 

r n~r 

all vertices w adjacent to v lie in the hyperplane xi + ■ ■ ■ + Xr = r — 1 since v — w = — 
where s G {1, . . . ,r} and I ^ {1, . . . , r}. □ 

In light of Theorem 62 we note that any triangulation of {v,Adj(v)} is unimodular 
where v E Iiic^Bm)- This follows from the fact that all adjacent vertices of v are in an 
n — 2 hyperplane, and thus any full- dimensional simplex X must contain v which implies 
G{X) is connected. 



67 



CHAPTER 4 

Applications to Optimization Through the Structure of 

Matroid Polytopes 

Here we present several heuristics and algorithms for non-linear matroid optimization 
and convex matroid maximization problems along with computational results. We give a 
description of our implementation of a speed-up of the algorithm proposed in [Onn03] for 
convex maximization and two primal heuristics, namely Local Search and Tabu Search, for 
convex minimization. We also give two heuristics for enumerating a subset of the projected 
bases. 



Before we start our description of the algorithms and heuristics, we remark that all of 
them rely on the geometry of the 1-skeleton graph of the matroid polytope [Sch03] . Recall 
Lemma 17 from Chapter 1, which states that the edges correspond to bases exchanges. 
These graphs have a lot of special structure. For example, these graphs are always Hamil- 
tonian [HH72] , and it is known that each two-dimensional face of every matroid polytope 
is either a triangle or a quadrilateral which means the graphs obtained when projecting are 
quite dense and easy to traverse [BCW07]. 



Given vi, . . . , v^, we define cone(vi, . . . , v^) := I Yli=i ^i^i I ^ f • If "P ^ R' 



is a polytope, we say two vertices vi,V2 G V are adjacent if they are contained in a 
one-dimensional face. Next we give a vital but elementary proposition necessary for our 
algorithms. 

Proposition 76. Let V OW^ be a polytope, vi, V2 vertices ofV, and W G M"'^". Then 
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V adjacent to vi 
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Proof. Let vi, . . . , be adjacent to vi. It follows from convexity of V that there 
exists Ai, . . . , Afc > such that 

V2 = Ai(vi - vi) H h Afc(vk - vi) + vi 

=^ Ww2 = AiVF(vi - vi) H h AfcM^(vk - vi) + Wwi. 

□ 

Recall from Chapter 1 the definitions WVm := { T^e^ | B e } C R'^, ij^WVjvi 
to mean the cardinality of WVm-, and WVm '■= conv(TyP;vj). Proposition 76 implies the 
following lemma. 

Lemma 77. Let V C M" be a polytope, v a vertex ofV, and W G M°'^". Then W^r e 
relint(VF'P) if and only if cone(Wv — Wv \ v adjacent to v). 

Define Adj(i?) := {B £ Bm I — = — ej }. Using Lemma 17 we get the following 
corollaries of Proposition 76 and Lemma 77. 

Corollary 78. Let M be a matroid on the ground set [n], W £ M*^^", B,B £ Bm- 
Weg = WeB+ ^ XB'iWeB' -Wbb) 

_B'eAdj(B) 

where Xb' > 0. 

Corollary 79. Let M be a matroid on the ground set [n] and W G M*^^" and B G Bm- 
Wbb e relint WV{M) if and only if cone ( Wbb' -WbbI B' e Adj(5) )=R'^. 

In order to prove that some of our heuristics are efficient, it is necessary to consider the 
generalized unary encoding of the weightings W G M'^^" [BLOW08]. We consider weights 
Wij of the form 

p 
k=l 

where p > 1 is a fixed integer, the distinct positive integers are binary encoded, and the 
integers 5ij are unary encoded. We can think of 6^ := G ^^x" for 1 < A; < p as 

unary-encoded matrices, with binary encoded a G where > for 1 < k < p. Then we 
have W := ((M^ij)) = X]fc=i (^k^^ ■ The generalized unary encoding includes three special 
cases: 



4.1. DESCRIPTION OF THE ALGORITHMS AND HEURISTICS 



69 



1. Unary-encoded weights: With p = 1 and oi = 1, we get the ordinary model of 
unary-encoded W = 6^. 

2. Binary-encoded {0, ai, . . . , ap}-valued weig hts: With 6'' > for ah k, and 
Yll:=i^ij ^ foi' ^iJi get the case of all Wtj in the set {0, ai, . . . , Op} 
having binary-encoded elements. 

3. 0\l-valued weights: This is, of course, the important common special case of 1 
and 2. 

With W encoded in this way it is possible to bound the size of #WVm- 

Lemma 80. Let d and p he fixed, M a matroid on n elements, and a generalized unary 
encoded matrix W := {{Wi,j)) = X]fc=i '^k^^ ^ Z"^^". Then #WVm is polynomially bounded 
in n, the binary encoding of ak, and the unary encoding of 6^ , for 1 < k < p. 

Proof. Let f3k = maxjj S'^j, for 1 < /c < p. Then 

WVm = { Wes \BgBm} 

^ < Xkak \\£{0,l,...,f3k rank(7W)}P, = rank(A^) \ 



k=l k=l 



□ 



The previous lemma bounds the size of #WVmi Y^t the number of integer points in 
the smallest rectangular region containing WVm may be exponential in the input. Some 
of our heuristics depend on the number of integers points in the smallest rectangular region 
containing WVm to be polynomially bounded in order to prove efficiency. 

Lemma 81 ([BLOW08]). Let d and p be fixed, M a matroid on n elements, and W G 
Z*^^". Then #WVm H Z'^ is polynomially bounded in n, and the unary encoding ofW, for 
l<k <p. 

Proof. Then 
WVm = { Wee \ B e Bm} 

{p p 
^Afc I AE{0,l,...,/3fcrank(A^)}P, ^Afc = rank(A^) I . 
k=l k=l J 
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□ 

4.1.1. Local and Tabu Search. First we present Heuristic 82 that starts at any base 
of Ai and proceeds by pivoting to its neighbors on V^M) as long as the pivot decreases the 
given function / through the weighting W. 

In the majority of our experiments / is concave or convex. In the case where / is 
concave, there will be a minimum on the boundary of WVm- If / is convex, then the 
minima may be in the interior of WVm ■ We emphasize that Heuristic 82 is not guaranteed 
to terminate at an optimum. 

Heuristic 82 (Local Search). 



hS{M,WJ,B) 

Input: Matroid M on n elements, B ^Bm,W ^ W^^"^, f -.W^ 
Output: A base B' G Bm such that f{We^) > f{WeB') G Adj(S') 

1: B' := B 

2: repeat 

3: OLDB' := B' 

4: if f{We^) < fiWes') for some B £ Adj(B') then 

5: B' := B 

6: until B' = OLDB' 

7: return B' 



If the objective function / is linear, then Heuristic 82 follows the non-degenerate steps 
of the simplex method. The following is an application of the well-known fact that the 
simplex method with a linear objective will terminate at a minimum [Sch86b]. 

Lemma 83. Let Ai be a matroid on [n], B G Bm, W £ M"^^", and /(x) = c • x, where 
c e R'^. LS(A^, W,g,B) terminates at B where g(Weg) = min{ g(WeB) \ B G Bm }• 

Lemma 83 is used by MOCHA in the implementation of our heuristics. For example, 
we can easily compute a tight rectangular region containing WV. We also use Lemma 83 to 
obtain an integer point on the boundary of WVm^ a prerequisite for some of our heuristics. 
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In section 4.3, we will show the practical limitations of Heuristic 82 and we give a 
theoretical bound on the running time of the previous algorithm given restrictions on W 
and d. 

Lemma 84. For every fixed d and p, Ai a matroid on [n] given by an independence 
oracle, a generalized unary encoded matrix W := ((VFjj)) = X]fc=i '^fc'^'^ ^ Z"'^", and f 
given by comparison oracle, then Heuristic 82 terminates in time polynomial in n, the 
binary encoding of ak, and the unary encoding of 6^, for 1 < k < p. 

Proof. For any B £ Bjn we can enumerate Adj(B) in polynomial time via an oracle 
for M by testing which of { B\{i} U {j} \ i £ B, j £ [n]\B } are bases. Thus, we can pivot 
in polynomial time in the input. Since Heuristic 82 always pivots to an adjacent base that 
is smaller, through / and W, than the current base, each point of WVm will be visited at 
most once. By Lemma 80, H^WVm ^ '^'^ is polynomially bounded in the input. Also, the 
arithmetic in Heuristic 82 is polynomial in n and the binary encoding of Ofc. □ 

Tabu search was first presented in [GI086] and has been widely used in combinatorial 
optimization. It begins at base B G Bm ^-nd pivots to an adjacent base B' if it is smaller than 
some other adjacent base, through / and the weighting W . This differs from Heuristic 82 
in that we allow pivots that are not necessarily smaller than the current base, through / 
and the weighting W. We record the smallest value /(VFe^) encountered and terminate 
after L pivots with no update to the minimum encountered. As for Heuristic 82, for the 
majority of our experiments / is concave or convex. We emphasize that Heuristic 85 is not 
guaranteed to terminate at the optimum for general /. 

Heuristic 85 (Tabu Search). 



TS(7W,5,VF,L) 

Input: Matroid M on [n], B £ Bm, W £ M°'^", L G N. 
Output: B' £ Bm 

1: B' := B 

2: VIS := {B} 

3: CURMIN := f{WeB'). 
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4: repeat 

5: if f{We^) < f{Weg) for some B,B e Adi{B')\ YIS then 

6: VIS := VIS U{B} 

7: B':=B 

8: if fiWes') < CURMIN then 

9: CURMIN := f(WeB'). 

10: until L pivots with CURMIN unchanged. 

11: return B' 



Heuristic 85 can also be modified such that we mark Wes as visited, not B, and only 
pivot to B such that Wes is unvisited. In this way, we could also prove that the running 
time of the modified algorithm is polynomial in the input, given the same restrictions for 
Lemma 87. Although, restricting to visiting new Wes could limit the possibilities of quickly 
converging to the optimum though. 

We observed through our computational experiments that Heuristic 82 and Heuristic 85 
work very well when optimizing concave or convex functions / over projected matroid 
polytopes as we will show in section 4.3. However, when / is arbitrary then the previous 
algorithms may not not reach the optimum. 

4.1.2. Listing Heuristics. In the remainder of the present section, we describe algo- 
rithms and heuristics aimed at tackling problems with an arbitrary balancing function / or 
using Pareto or minmax optimization. The following algorithms do not explicitly evaluate 
/ but are presented as listing algorithms. Note that if we have all of the projected matroid 
bases then it is simple to extract the Pareto optima through a straightforward pairwise 
comparison. Likewise for optimizing over / or with a min-max objective, we simply evalu- 
ate over all projected bases found. This is a similar methodology as [BLOW08] where the 
authors prove efficient deterministic algorithms, given sufficient conditions on the input, to 
list all projected bases. 

In [BLOW08] the authors used matroid intersections to solve the problem: Given x G M*^ 
find B £ Bm such that Wes = x if such a B exists. Guided by the success of our previous 
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heuristics, we had the new idea of using either Heuristic 85 or Heuristic 82 to find a i? E Bm 
such that Wes = x. The novelty of our heuristic is that we use fast convex optimization 
as a subroutine to solve no n- linear and Pareto optimization problems. 

Given an x' € Z*^, Heuristic 86 forms a convex function /x' that minimizes at x' and 
/x'(x') = 0. Our heuristic calls Heuristic 82 or Heuristic 85 with the function /x' and if it 
returns a base B £ Bm such that /x'CVFe^) = then Wes = x. 

Heuristic 86 (Pivot Test). 



PT{M,B,W,t,S) 

Input: Matroid M on n elements, B £ Bm, W £ M'^''", t G N, finite set S C Z"*. 
Output: PT C Bm such that {Wbb \ B £ VJ} O S. 
1: PT := 

2: for each x' € S" do 

3: /.'(x):=Etl(^-^i)' 

4: for 1, . . . , t do 

5: B' := random basis of M 

6: B := LS(7W, ^, fx',B') (Or use Heuristic 85) 

7: if /(lyes) = then 

8: PT := PT U S 

9: Break for loop 

10: return PT 



Lemma 87. For every fixed d and p, Ai a matroid on [n] given by an independence 
oracle, a generalized unary encoded matrix W := {{Wij)) = X]fc=i o-k^^ £ Z*^^", a finite set 
S C Z'^, then Heuristic 86 terminates in time polynomial in n, t, the size of S, the binary 
encoding of ak, and the unary encoding of 6^, for 1 < k < p. 



Proof. This follows from the fact that Heuristic 82 and Heuristic 85 are polynomial 
under these assumptions. Moreover, we call Heuristic 82 or Heuristic 85 at most t\S\ 
times. □ 
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Typically we wish to use Heuristic 86 to try to enumerate a subset of W-p. Thus, if we 
further assume W is unary encoded, then the smallest box containing WVm is polynomial 
which we use as our input S in Heuristic 86. 

For general non-linear multi-criteria matroid optimization problems, when the function 
/ we are minimizing is concave, then some optimum will be a vertex of WVm- Thus, 
in such a case, it is sufficient to enumerate the vertices of WVm t^st / over those 
points. Furthermore, some vertices are also Pareto optimal, but in general, the vertices 
are not all the Pareto optima. Later we will give an algorithm that uses the vertices of 
WVm to facilitate finding other Pareto optima. We will prove that the following algorithm 
will enumerate all vertices of WVm possibly other integral points on the boundary of 
WVm- Algorithm 88 starts at a base B £ Bm such that We^ is on the boundary of WVm- 
If B' € Adj(i?), Wes' is a newly seen point and is on the boundary of WVm^ record B' 
and We^', pivot to B', and continue. 

Algorithm 88 (Projected Boundary). 
FB{M,W,B) 

Input: Matroid A4 on n elements, W G M'^^", B G Bm such that Wes is an extreme 
point of WVm - 

Output: CH C Bm such that {Wes \ B £ CE} ^ vert(WT^). 
1: CH := {B'} 

2: Mark B' as unvisited in CH 
3: PB := {Wes'} 
4: repeat 

5: B' := first unvisited base in CH 
6: Mark B' as visited. 
7: for B G Adj(5') do 

8: if Weg is an extreme point of WVm then 
9: if Weg ^ PB then 

10: CH := CH U {B} 
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11: Mark B as unvisited. 

12: PB := PB U{M^e^} 

13: until PB unchanged. 
14: return CH 



In section 4.3 we will see that Algorithm 88 works well. We believe this is due to the 
following: the projected bases are often highly clustered and with much fewer projected 
bases near the boundary of WVm- To see how large l^~^x, for x G M'^, can be, see Figure 
1 in [GLM08] . The following lemma proves that the output of Algorithm 88 will contain all 
the vertices of the convex hull of WVm ■ 

Lemma 89. Let M. be a matroid on [n], W G M^^", B e Bm such that Wbb IS an 
extremal point of WVm ■ Then 



Yeit{WVM) ^ W{FB{M, W,B))C {WVm n OWVm)- 

Proof. Let B',B G Bm where Wes' is an extremal point and Weg is a vertex of 
WVm ■ Moreover let We b' and We^ be on a 1-face F of WVm ■ By Corollary 78 there exists 
B" £ Adj(i?') such that either We^" = We^ or We^" is a positive convex combination 
of Wgb' and We^, i.e. Wgb" G F. Thus Wgb" is extremal and Algorithm 88 can always 
pivot towards every vertex of WVm a-iid will output all vertices of WVm ■ ^ 

B 

\ 

Figure 4.1. A pivot from B' to B" towards B in Algorithm 88. 

In [Onn03] it was shown that the number of vertices of WVm is polynomially bounded 
when the number of weightings d is fixed (see Proposition 2.1 and Lemma 2.3). Moreover 
the author showed that the vertices of WVm can be found in polynomial time. A similar 
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result can be found in [BLOW08]. Algorithm 88 can (and does in practice) pick up non- 
vertex integral boundary points. Hence, a bound on the number of vertices is not sufficient 
to guarantee efficiency of the algorithm. 

Lemma 90. For fixed d and p, M a matroid on [n] given by an independence oracle, 
a generalized unary encoded matrix W := {{Wij)) = Ylk=i^k^^ ^ Z'^^"', Algorithm 88 
runs in time polynomial in n, the binary encoding of au, and the unary encoding of 5^ , for 
l<k<p.. 

Proof. Since #WVm n Z'^ is polynomially bounded and 7^ Adj(i?) < we only have 
to show that deciding if We^ is an extremal point of WVm is polynomial in the input. 
This is equivalent to deciding if cone{WeB — We^ \ B £ Adj(i?) } is pointed. We can 
decide this by solving the feasibility problem: 

|x,-x G cone{WeB -We^ \ B G Adj(i?)} | x > o} 

□ 

We will show in section 4.3 the effectiveness of Heuristic 86 (PT) in enumerating all 
projected bases. When finding Pareto optima though it is not necessary to enumerate all 
projected bases. The following heuristic combines Algorithm 88 and Heuristic 86 to find 
a set of points which could be Pareto optima of WVm- We observe that the min-max 
optimum is also a Pareto optima. Hence Algorithm 88 and Heuristic 86 could return the 
min-max optimum of WVm ■ The main idea is that if we have the boundary points, then the 
regions where other Pareto optima lay is determined by the convex hulls of points derived 
from a triangulation of the boundary. 

Heuristic 91 (Boundary and Triangular Region Pareto Test). 



BTRPT(7W,VF,t) 

Input: Matroid M on n elements, W G M''^", t G N. 
Output: PO C Bm 

1: PO := PB{M,W) 

2: for i3 G PO do 
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3: for B' ePO\B do 

4: if WeB' >Pareto Wsb then 

5: PO := PD \ B' 

6: PPO := WPO projected on the hyperplane H := {x \ J2^i — 

7: Triangulate PPO. 

8: for each facet F of PPO do 

9: TRIPOINTS := {Bases in PO corresponding to vertices of F} 
10: for i = 1, . . . d do 

11: X™'"^' := max( (WeB)i \ B £ TRIPOINTS ) 

12: PO := PT(>i,S, W,t,conv(TRIPOINTS U X''^'^'')) U PO 

13: for S G PO do 

14: for B' e PO\B do 

15: if Wbb' >Pareto Wbb then 

16: PO := PO \ B' 

17: return PO 




Except for the boundary points, Heuristic 91 is not guaranteed to return points that 
are Pareto optima since Heuristic 86 (PT) may miss some projected bases in the test regions. 



4.1. DESCRIPTION OF THE ALGORITHMS AND HEURISTICS 



78 



Motivated by the need of optimizing any function / or exploring all Pareto or minmax 
optima, we devised a variation of the breadth-first search algorithm that would limit the 
search to some depth of the tree. Our intuition was that the graph of base exchanges is 
highly connected, thus small depth was all that is necessary for listing a large proportion 
of all projected bases, as we will see in section 4.3.: 

Heuristic 92 (Different Fiber BFS). 



BFBFS{M,W,B,d,l,PB) 

Input: Matroid M on n elements, W G M'^'<", B £ Bm, d, I £ N, PB <Z WVm- 
Output: DFBFS ^ {Wes \ B £ Bm} 

1: if l>d then 

2: return 

3: PB := PR U{WeB} 

4: PTT := 

5: for each B' G Adj(S) do 
6: if WeB' i PB then 
7: if Wgb' / Wgb then 
8: PB := PB U VFe^/ 

9: PTT := PTT U B' 

10: for each B' G PTT do 

11: PB := PB U DFBFS(7W,W^,B',d,/ + 1,PB). 
12: return PB 



Heuristic 92 begins at a base B G Bj^. It adds the projected point We.B to the set PB, 
short for projected bases. It then enumerates all adjacent bases Adj(i3) of B. If a neighbor 
i3''s projected point Wbb' ^ PB and WeB ^ WeB', then we add Wes' to PB and recursively 
call DFBFS on B' . We allow a parameter d which determines the recursive depth allowed. 
Heuristic 92 takes its name, different fiber breadth-first-search, from the fact that we do 
not allow pivots that evaluate the same under the weighting W. As a small consequence. 
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this guarantees that the number of times DFBFS is called is bounded by H^WVm ^ useful 
fact if the input is bounded as in Lemma 80. 

4.2. Software Implementation 

Our heuristics are implemented in C++ and take advantage of the object-oriented 
paradigm. The software MOCHA^[DHLO09] reads in either a vectorial matroid, represented 
by an m X n floating-point matrix or a graphical matroid, represented by an n x n adjacency 
matrix. The weightings are read in as a d x n floating point matrix. 

For vectorial matroids, the rank of a subset of columns is computed using LAPACK, a 
standard and robust linear algebra package [ABB+99]. For vectorial matroids with elements 
in Z, MOCHA has the option of using the GMP arbitrary-precision software to perform 
Gaussian elimination using exact arithmetic [Gea09]. Enumeration of the neighbors of a 
base B is done by calculating the rank of B\i U j for all i £ B and j € If it is full 

rank then it is returned as a neighboring base. Random bases are determined by randomly 
choosing a rank(A^) sized subset A C [n] and checking if rank(74) = rank(7W), repeating 
until such an A is found. 

For graphical matroids, the rank of the matroid, and any A C [n], is determined by 
calculating the size of a spanning forest via breadth-first-search. To enumerate the neighbors 
of G Bm we first calculate all paths in B using dynamic programming. Adding any 
element (edge) j ^ B to B will create a cycle C. We use our pre-calculated all paths of B 
to quickly determine C. Then all subsets B\i U j where i £ C will be an adjacent base to 
B. This may not be the optimal solution to graphical matroid adjacency enumeration, but 
it is straight-forward. 

The Projected Boundary algorithm (Algorithm 88) is only implemented for d = 2. This 
is due to the computational expense of determining if We^ is an extreme point of WVm 
on line 8. When d = 2 this is easy to check by sorting the projected rays between B and its 
neighbors. If there exists two sorted rays with angle larger than vr then We^ is extreme. 

For correctness, we compare the number of projected spanning trees found using our 
methods versus the real total number of projected spanning trees. We used an algorithm for 
generating all of the spanning trees in undirected graphs presented by Matsui [Mat97] . The 

-'^Matroid Optimization Combinatorial Heuristics and Algorithms 



4.3. COMPUTATIONAL RESULTS 



80 



algorithm requires 0{n + m + rn) time when the given graph has n vertices, m edges, and 
r spanning trees. For outputting all of the spanning trees explicitly, this time complexity is 
optimal. We also implemented the asymptotic 0\1 polytope vertex-estimation presented in 
[BS07a]. This gives us the ability to estimate the number of bases of matroid polytopes in 
order to better understand the proportion of bases to projected bases for problems where 
full enumeration is intractable. 

4.3. Computational Results 

Now we present our experiments for which we used six roughly comparable machines 
(see Table 4.38). For our tests. Heuristic 82 uses the pivot rule that fiWe^) < /(VFe^/) 
for all B G Adj(-B') in line 4. Heuristic 85 uses the pivot rule that /{WBj^) < f{Weg) for 
all B,B e Adi{B') \ VIS in line 5. 

4.3.1. Calibration Set. Our first goal was to perform experiments on matroids for 
which we can compute all bases in order to better understand our heuristics and algo- 
rithms. We generated fifteen connected random graphs: gn9el8, gn9e27, gnl0e22, gnl0e28, 
gnlOeSS, gnllelS, gnlle20, gnlle27, gnlle41, gnl2el6, gnl2e24, gnl2e33, gnl3el9, gnl3e29, 
gnl3e39 which we will refer to as our calibration set. The names of our graphs follow the 
simple nomenclature 5n[7^ nodes] e[7^edges] . We consider two, three and five criteria, i.e. 
number of weightings. We further consider three different ranges of integral weights for 
each criteria. For the calibration set we adopt the following nomenclature 

5n[#nodes] e[#edges] d[7^criteria] w;[low weight] u;[high weight] 

where we generated random integral weightings between [low weight] and [high weight] . 
First we simply compare the number of spanning trees of our calibration set to the number 
of projected spanning trees. 

• Table 4.1 shows the calibration set with two weightings (criteria) and integral 
weights - 20, - 100, and - 1000. 

• Table 4.2 shows the calibration set with three weightings (criteria) and integral 
weights - 20, - 100, and - 1000. 

• Table 4.3 contains the calibration set with five weightings (criteria) and integral 
weights — 1,0 — 2, and — 5. 
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We give the exact number of projected spanning trees and compare versus the exact number 
of spanning trees. 

4.3.1.1. Calibration Set - Heuristic 92. There are four parameters to our implementa- 
tion of Heuristic 92 (DFBFS): number of searches A^, BFS depth, boundary retry limit 
and random retry limit. First we attempt to find a new boundary projected bases using 
Heuristic 82 (Local Search) and a random direction. For every new projected base, we run 
Heuristic 92 (DFBFS) with the given depth parameter. We attempt to find N new bound- 
ary projected bases and give up if we exceed the boundary retry limit. Next our algorithm 
will generate a random base and project it by the weighting. If it is a new projected base 
we run Heuristic 92 (DFBFS) with the given depth parameter. We attempt to find N new 
random projected bases and give up if we exceed the random retry limit. 

• Table 4.4 shows Heuristic 92 (DFBFS) on the calibration set in two criteria. 

• Table 4.5 shows Heuristic 92 (DFBFS) on the calibration set in three criteria. 

• Table 4.6 shows Heuristic 92 (DFBFS) on the calibration set in five criteria. 

We list the machine used, seconds required, number of searches, BFS depth, boundary retry 
limit and the random retry limit. We noticed a decrease in the proportion of projected trees 
found to all projected trees as the dimension and weights increased. It is also important 
to note, as seen in Table 4.1 and Table 4.2, as the dimension and weights increase, the 
proportion of all projected trees to all spanning trees increases. 

In our experiments, in almost all cases. Heuristic 92 (DFBFS) terminated by exceed- 
ing the boundary and random retry limit and not the number of searches. This can be 
attributed to the phenomenon that bases are projected in a relatively tight area, with few 
projected bases near the boundary. Other conditions for picking the initial base for Heuris- 
tic 92 (DFBFS) could yield better results. For example, keeping a random leaf of the 
previous truncated BFS. Or we could compute the boundary of the projected bases and 
attempt to find a new initial base using Heuristic 86 (PT) on the "holes" in the convex hull. 

4.3.1.2. Calibration Set - Heuristic 86. Our implementation of Heuristic 86 (PT) uses 
Heuristic 82 (LS) with the 2d directions { itej | i £ {1, . . . ,d}} to compute a bounding box 
containing all projected base. 
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Figure 4.3. 96.18% projected bases of gnlle41d2wOwlOO using Heuris- 
tic 92 (DFBFS) 



• Table 4.7 shows the result of Heuristic 86 (PT) on our calibration test set with two 
criteria and weights — 20 and — 100 using Heuristic 82 (LS) as our test method 
subroutine. We give the proportion of projected trees found versus all projected 
trees. Naturally with larger weight values, the projected trees will be contained in 
a larger box, requiring more time. 

• Table 4.8 shows the result of Heuristic 86 (PT) on our calibration test set with 
two criteria and weights — 20 using Heuristic 85 (TS) as our test method sub- 
routine. The run times using Heuristic 85 (TS) are longer than Heuristic 82 (LS), 
but we find nearly all the projected trees. The worse case is the sparse graph 
gnllel3d2w0'w20.mo where we find only 90.24% of all projected trees. 

Heuristic 86 (PT) has two major advantages; 
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• Heuristic 86 uses very little memory since Heuristic 85 (TS) and Heuristic 82 (LS) 
use little memory; 

• Heuristic 86 can be distributed since one can partition up test regions and run 
each region as a separate instance. 

4.3.1.3. Calibration Set - Algorithm 88. Table 4.9 shows the number of projected span- 
ning trees found on the boundary using Algorithm 88 (PB). All computations took less 
than one second. To find a starting base, Algorithm 88 (PB) generates a random direction 
and calls Heuristic 82 (LS). It is interesting to note that the graph gnl3e39d2w0'wl000 has 
15,037,589 projected spanning trees, yet there are only 26 points found by Algorithm 88 
(PB). 

4.3.1.4. Calibration Set - Heuristic 91. Table 4.10 shows the results of Heuristic 91 
(BTRPT) on the calibration set with weightings — 20 and — 100 in two criteria. We use 
Heuristic 85 (TS) as our subroutine in Heuristic 86 (PT). We give the seconds, the exact 
number of Pareto optima and the number of Pareto optima found. Heuristic 91 (BTRPT) 
not only found the correct number of Pareto optima, it found the correct Pareto optima in 
all cases. Heuristic 91 (BTRPT) can be distributed: once the boundary is computed and 
the regions are found. Heuristic 86 (PT) can be run as separate instances for each region 
(or subdivided further). 

4.3.1.5. Calibration - Local and Tabu Search For Convex Minimization. For our exper- 
iments of Heuristic 82 (LS) and Heuristic 85 (TS) on the calibration set we minimized over 
two convex functions: (x — x)'^ and (x — x)^. First we choose x G WVm such that it is 
an interior point of WVm- Second we consider x to be rational non- integer point which 
is an interior point of WVm- For the integral case, we can easily detect if we are at the 
minimum since our objective will evaluate to zero if so. For the rational case, we verify the 
global minimum by evaluating our objective on all projected spanning trees. For all tests 
we perform 1000 minimizations with random starts and record the number of successes. 
For easy reference, in the cases of Heuristic 82 (LS) we also repeat the exact number of 
spanning trees and exact number of projected spanning trees and the proportion (given as 
percentage) between them. 
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Figure 4.4. Green Pareto optima and blue test regions of 
gnlle41d2w0wl00 using Heuristic 91 (BTRPT) 



• Table 4.12, Table 4.13, Table 4.14 show Heuristic 82 (LS) where the minimum is 
the integer point described above minimizing (x — x)^. 

• Table 4.15, Table 4.16, Table 4.17 show Heuristic 82 (LS) where the minimum is 
the integer point described above minimizing (x — x)^. 

• Table 4.24, Table 4.25, Table 4.26 show Heuristic 82 (LS) where the minimum is 
the rational point described above minimizing (x — x)'^. 

• Table 4.27, Table 4.28, Table 4.29 show Heuristic 82 (LS) where the minimum is 
the rational point described above minimizing (x — x)^. 

• Table 4.18, Table 4.19, Table 4.20 show Heuristic 85 (TS) where the minimum is 
the integer point described above minimizing (x — x)'^. 

• Table 4.21, Table 4.22, Table 4.23 show Heuristic 85 (TS) where the minimum is 
the integer point described above minimizing (x — x)^. 

• Table 4.30, Table 4.31, Table 4.32 show Heuristic 85 (TS) where the minimum is 
the rational point described above minimizing (x — x)^. 
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• Table 4.33, Table 4.34, Table 4.35 show Heuristic 85 (TS) where the minimum is 
the rational point described above minimizing (x — x)'*. 

We first observe that the success of Heuristic 82 (LS) and Heuristic 85 (TS) decrease 
as the proportion of projected spanning trees to all spanning trees increases. Second, for 
Heuristic 85 (TS) we see a noticeable increase in the number of successes as the Tabu Limit 
ranges from 1 to 10 indicating that in many cases, a low Tabu Limit is sufficient to reach 
the global minimum. 

4.3.2. Pushing the Limits. 

4.3.2.1. Sparse Graphs - Solids. We now compare sparse graphs versus the dense graphs 
presented in the calibration set. We chose fifteen planar graphs of the 1-skeleton of three- 
dimensional polytopes. Table 4.11 shows the result of Heuristic 92 (DFBFS) on these 
examples. For some of these graphs, the number of spanning trees can be explicitly enu- 
merated, hence we give the exact number of projected spanning trees. When the weights 
take value from — 1000 we omit the last five graphs as DFBFS exceeded over a day of 
computation or the machines (/uz2;j/, truth) ran out of memory. 

4.3.2.2. Experimental Design. In [BLMA"'"08] the authors first proposed using non- 
linear matroid optimization to solve the statistical experimental design problem, which 
is explained in detail in [FH80] . The experimental design problem can be briefly described 
as attempting to learn an unknown system whose output y is an unknown function <I> with 
input x = (xi, . . . ,Xk) G M'^. To learn the system, experiments are performed using input 
Pi = {pi^i, . . . ,Pi,k) and the output = $(pi) is measured. Then, based on the experiments 
we wish to fit a model for the system, that is, determine an estimation <I> of the function $ 
such that it: lies in a prescribed class of functions, is consistent with the outcomes of the 
experiments, and minimizes the aberration (some suitable criteria) among models in the 
class. 

A deterministic polynomial time algorithm was proved in [BLMA+08] which solves non- 
linear matroid optimization over arbitrary matroids (presented by evaluation oracle) when 
the number of weightings are fixed and the weightings are a general unary encoding. If d is 
the number of weightings, the weights only take value on p many integers, and the matroid 
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is on n elements, then their algorithm for arbitrary matroids uses the matroid intersection 
algorithm nP'' a polynomial number of times. 

However, for vector matroids they provide a faster algorithm (compared to their algo- 
rithm for arbitrary matroids) for unary encoded weights. If uj is the size of the maximum 
unary encoded element of the weights, then this improved algorithm requires solving a 
(rank(A^)u; + l)'^ x (rank(A4)6<-' + l)'^ Vandemonde equation. This non-linear vector matroid 
optimization algorithm has been implemented on IBM's Blue Gene/L computer [GLM08]. 
At the authors disposal was the 360TF, 128K processor Blue Gene/L supercomputer at 
Lawrence Livermore National Laboratories. We present now two instances P20 and P28 
which were solved using the deterministic vector matroid algorithm on said supercomputer. 
Each is an instance of an experimental design problem encoded as a vector matroid non- 
linear optimization program. Both their algorithm and our heuristics (Heuristic 92 and 
Heuristic 86) do not optimize a particular objective function, but instead list the projected 
bases (note our heuristic is not guaranteed to find all projected bases). 

Table 4.36 shows the size of the instances P20 and P2S, the true number of projected 
bases and the number of seconds required by Blue Gene/L to find the projected bases. We 
also show two runs of Heuristic 92 (DFBFS) using Fuzzy with the number of projected 
bases found and seconds required. What is noteworthy is that Heuristic 92 found all the 
projected bases in at least one of the two runs. Moreover, the other runs of Heuristic 92 
found nearly all the projected bases. This demonstrates that our heuristics are useful when 
solving these difficult problems, especially in light of the fact that they were computed using 
a relatively cheap computer system. An additional point is that the deterministic vector 
matroid algorithm requires very high precision in order to solve the large Vandemonde 
systems derived from the matroid optimization problem. On the other hand, our heuristics 
only depend numerically on solving at most rank(A4) x rank(7W) linear systems. In fact, our 
software has an option to use exact arithmetic, albeit with a slow down in solution times. 
We emphasize the drastic difference in running times and computational power between 
Blue Gene/L and Fuzzy. 

In Table 4.37 we show the result of Heuristic 86 (PT) on i-20- We divided up the 
region of feasible projected bases into twelve disjoint regions and ran twelve instances of 
Heuristic 86 (PT) concurrently. We report the times and the number of projected bases 
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found for each region. We also give the total number of projected bases found which is the 
exact number of projected bases. For the search in Heuristic 86 (PT) we used Heuristic 85 
(TS) 10 times on each point with a retry limit of 20. We performed this experiment to 
exhibit the effectiveness and ease by which Heuristic 86 can be distributed. 

4.3.3. Discussion. Non-linear matroid optimization problems are very difficult. They 
are NP-complete in general and exponential in some cases [BLMA"^08]. Our goal was to 
present and explore the practicality of new heuristics and algorithms for solving these prob- 
lems. The effectiveness of our new techniques rely on two important properties: traversing 
using the local adjacency inherent in matroids is easy, and that although there many be 
exponentially many bases, the number of projected bases is manageable. 

The fact that many of our heuristics are not guaranteed to find the optimal solution or 
list all the projected bases is countered by the fact that they are not overly complex and 
very fast. For instance. Heuristic 92 (DFBFS) is not guaranteed to find all the projected 
bases but in many of the tests it found a large portion of them. In practice, as seen in 
subsubsection 4.3.2.2, since Heuristic 92 (DFBFS) finds a large portion of the projected 
bases, we would simply run the heuristic many times and union the results. 

Most all of our heuristics and algorithms complement each other. For example, we can 
start by finding the vertices using Algorithm 88 (PB) then enumerate as many projected 
bases as possible using Heuristic 92 (DFBFS). Next we can run Heuristic 86 (PT) on any 
regions we suspect there to be more projected bases. 

Although we do not prove it theoretically, our heuristics and algorithms use minimal 
memory. For instance, the driving part of all of our heuristics and algorithms is enumerating 
neighboring bases, which is quadratic in the number of elements of the matroid. The biggest 
memory hog is Heuristic 92 (DFBFS) which at worst stores all the projected bases found 
thus far. 

This leads to the ease by which Heuristic 86 (PT) can be distributed. In fact, the 
only overhead would be combining the solutions, which would be a simple union. Ideally 
if we had X machines with Y processors, we can first find a box containing the potential 
projected bases using 2n linear programs using Heuristic 82 (LS). Next we partition the 
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region of potential projected bases into XY regions and run an instance of Heuristic 86 
(PT) for each. 

Table 4.1. Calibration test set: Exact #Spanning trees vs. ^projected spanning trees in 2 criteria. 
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Table 4.2. Calibration test set: Exact #Spanning trees vs. #projected spanning trees in 3 criteria. 
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Table 4.3. Calibration test set: Exact #Spanning trees vs. #projected spanning trees in 5 criteria. 
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Table 4.7. Calibration test set: Heuristic 86 (Pivot Test) using Heuristic 82 (Local Search), 
10 searches per point. 



B 


odes 


tu 
bo 
"0 


imension 




^eights 


omputcr 




sconds 




Projected 
Trees 




Trees Found 
by DFBFS 


ercent 


2; 


2; 


H 


Q 






o 












* 




gn9cl8d2w0w20 


9 


18 


2 


- 


-20 


mocrOl 




92 


1 


742 


1 


121 


64.35 




9 


27 


2 


- 


-20 


mocrOl 




169 


5 


346 


4 


939 


77.83 


gnl0G22d2w0w20 


10 


22 


2 


- 


-20 


mocrOl 




101 


3 


957 


3 


318 


83.85 




10 


28 


2 


- 


-20 


inocrOl 




207 


7 


131 


g 


163 


86.43 


gnl0c33d2w0w20 


10 


33 


2 


0- 


-20 


inocrOl 




356 


10 


833 


8 


842 


81.62 


gnllel3d2w0w20 


11 


13 


2 


0- 


-20 


inocrOl 




5 




41 




37 


90.24 


gnlle20d2w0w20 


11 


20 


2 


0- 


-20 


inocrOl 




89 


2 


173 


1 


328 


61.11 


gnllo27d2w0w20 


11 


27 


2 


0- 


-20 


inocrOl 




273 


7 


515 


5 


803 


77.22 


gnlle41d2w0w20 


11 


41 


2 


0^ 


-20 


inocrOl 




806 


16 


457 


14 


747 


89.61 


gnl2el6d2w0w20 


12 


16 


2 


0- 


-20 


inocrOl 




13 




154 




144 


93.51 


gnl2e24d2w0w20 


12 


24 


2 


0- 


-20 


mocrOl 




233 


5 


647 


3 


701 


65.54 


gnl2e33d2w0w20 


12 


33 


2 


0- 


-20 


inocrOl 




357 


9 


364 


7 


599 


81.15 


gnl3el9d2w0w20 


13 


19 


2 


0- 


-20 


inocrOl 




68 


1 


608 


1 


051 


65.36 


gnl3e29d2w0w20 


13 


29 


2 


0- 


-20 


inocrOl 




359 


8 


474 


6 


407 


75.61 


gnl3c39d2w0w20 


13 


39 


2 


0- 


-20 


inocrOl 


1 


160 


18 


468 


15 


082 


81.67 


gllcJdOU^WUWiUU 


9 


18 


2 


— 


100 


inocr02 




944 


2 


858 


\ 


346 


47.10 


gn9e27d2w0wl00 


9 


27 


2 


— 


100 


inocr02 


5 


923 


89 


092 


8 


760 


09.83 




10 


22 


2 


— 


100 


inocr02 


4 


315 


37 


204 


5 


249 


14.11 




10 


28 


2 


— 


100 


inocr02 


g 


224 


101 


334 


8 


236 


08.13 


gnl0c33d2w0wl00 


10 


33 


2 


0- 


100 


inocr02 


11 


925 


166 


427 


19 


247 


11.56 


gnllcl3d2w0wl00 


11 


13 


2 


0- 


100 


inocr02 




215 




41 




41 


100.00 


gnllc20d2w0wl00 


11 


20 


2 


0- 


100 


inocr02 


2 


737 


6 


580 


2 


107 


32.02 


gnlle27d2w0wl00 


11 


27 


2 


0- 


100 


mocr02 


6 


327 


81 


803 


6 


790 


08.30 


gnlle41d2w0wl00 


11 


41 


2 


0- 


100 


inocr02 


26 


144 


309 


961 


29 


373 


09.48 


gnl2el6d2w0wl00 


12 


16 


2 


0- 


100 


inocr02 




272 




162 




155 


95.68 


gnl2e24d2w0wl00 


12 


24 


2 


0- 


100 


inocr02 


11 


812 


92 


813 


11 


596 


12.49 


gnl2e33d2w0wl00 


12 


33 


2 


0- 


100 


inocr02 


15 


676 


192 


122 


16 


705 


08.70 


gnl3el9d2w0wl00 


13 


19 


2 


0- 


100 


inocr02 


1 


859 


3 


255 


2 


010 


61.75 


gnl3o29d2w0wl00 


13 


29 


2 


0- 


100 


inocr02 


15 


365 


164 


617 


12 


986 


07.89 


gnl3o39d2w0wl00 


13 


39 


2 


- 


100 


inocr02 


32 


958 


315 


881 


38 


530 


12.20 




9 


18 


2 


— 


100 


inocrOl 


3 


395 


2 


905 


1 


659 


57.11 


o-nQn97rl '^nrOii790 


9 


27 


2 


— 


100 


inocrOl 


27 


537 


134 


039 


16 


842 


12.57 


(rn 1 0n99r1 ^wrnw90 


10 


22 


2 


— 


100 


inocrOl 


18 


141 


42 


887 


8 


564 


19.97 




10 


28 


2 


— 


100 


inocrOl 


50 


429 


238 


529 


39 


819 


16.69 


gnl0c33d3w0w20 


10 


33 


2 


0- 


100 


inocrOl 


76 


051 


411 


730 


47 


606 


11.56 


gnllel3d3w0w20 


11 


13 


2 


0- 


100 


inocrOl 




271 




41 




40 


97.56 


gnlle20d3w0w20 


11 


20 


2 


0- 


100 


inocrOl 


6 


726 


6 


603 


2 


851 


43.18 


gnlle27d3w0w20 


11 


27 


2 


0- 


100 


inocrOl 


39 


970 


146 


672 


21 


790 


14.86 


gnlle41d3w0w20 


11 


41 


2 


0- 


100 


inocrOl 


216 


987 


959 


469 


120 


142 


12.52 


gnl2el6d3w0w20 


12 


16 


2 


0- 


100 


inocrOl 


1 


329 




162 




155 


95.68 


gnl2e24d3w0w20 


12 


24 


2 


0- 


100 


inocrOl 


32 


498 


111 


201 


16 


698 


15.02 


gnl2e33d3w0w20 


12 


33 


2 


0- 


100 


inocrOl 


107 


663 


470 


609 


53 


432 


11.35 


gnl3cl9d3w0w20 


13 


19 


2 


0- 


100 


inocrOl 


9 


980 


3 


358 


2 


255 


67.15 


gnl3c29d3w0w20 


13 


29 


2 


0- 


100 


mocrOl 


62 


506 


264 


949 


34 


109 


12.87 


gnl3e39d3w0w20 


13 


39 


2 


0- 


100 


mocrOl 


232 


975 


930 


322 


126 


044 


13.55 
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Table 4.8. Calibration test set: Heuristic 86 (Pivot Test) using Heuristic 85 (Tabu Search), 
10 searches per point, 20 search limit. 



Name 


Nodes 


Edges 


Dimension 




Weights 


Computer 




Seconds 




^Projected 
Trees 




#TrGcs Found 
by DFBFS 


Percent 


gn9el8d2w0w20 


9 


18 


2 


- 


-20 


mocr04 




333 


1 


742 


1 


643 


94.32 


gn9o27d2w0w20 


9 


27 


2 


- 


- 20 


mocr04 




518 


6 


346 


6 


299 


99.26 


gnl0c22d2w0w20 


10 


22 


2 


- 


- 20 


mocr04 




394 


3 


957 


3, 


878 


98.00 


gnl0c28d2w0w20 


10 


28 


2 


- 


- 20 


mocr04 




712 


7 


131 


7 


085 


99.35 


gnl0c33d2w0w20 


10 


33 


2 


- 


- 20 


mocr04 


1,146 


10 


833 


10 


739 


99.13 


gnllcl3d2w0w20 


11 


13 


2 


- 


- 20 


mocr04 




41 




41 




37 


90.24 


gnllc20d2w0w20 


11 


20 


2 


0- 


- 20 


mocr04 




440 


2 


173 


1 


979 


91.07 


gnllc27d2w0w20 


11 


27 


2 


0- 


- 20 


mocr04 




971 


7 


515 


7 


441 


99.02 


gnlle41d2w0w20 


11 


41 


2 


0- 


- 20 


mocr04 


2,484 


16 


457 


16 


379 


99.53 


gnl2el6d2w0w20 


12 


16 


2 


0- 


- 20 


mocr04 




85 




154 




153 


99.35 


gnl2c24d2w0w20 


12 


24 


2 


0- 


- 20 


mocr04 




898 


5 


647 


5 


449 


96.49 


gnl2c33d2w0w20 


12 


33 


2 


0- 


- 20 


mocr04 


1, 


082 


9 


364 


9 


270 


99.00 


gnl3cl9d2w0w20 


13 


19 


2 


0- 


- 20 


mocr04 




340 


1 


608 


1 


521 


94.59 


gnl3c29d2w0w20 


13 


29 


2 


0- 


- 20 


mocr04 


1, 


102 


8 


474 


8 


408 


99.22 


gnl3c39d2w0w20 


13 


39 


2 


0- 


- 20 


mocr04 


3, 


233 


18 


468 


18 


316 


99.18 


gn9el8d2w0wl00 


9 


18 


2 


— 


100 


mocr04 


6, 


231 


2 


858 


2 


556 


94.32 


gn9e27d2w0wl00 


9 


27 


2 


— 


100 


mocr04 


30, 


162 


89 


092 


31 


929 


99.26 


gnl0e22d2w0wl00 


10 


22 


2 


— 


100 


mocr04 


25, 


007 


37 


204 


20 


976 


98.00 


gnl0o28d2w0wl00 


10 


28 


2 


— 


100 


mocr04 


33, 


342 


101 


334 


50 


722 


99.35 


gnl0o33d2w0wl00 


10 


33 


2 


- 


100 


mocr04 


56, 


517 


166 


427 


82 


275 


99.13 


gnllel3d2w0wl00 


11 


13 


2 


- 


100 


mocr04 


1, 


561 




41 




41 


90.24 


gnlle20d2w0wl00 


11 


20 


2 


0- 


100 


mocr04 


17, 


385 


6 


580 


4 


972 


91.07 


gnlle27d2w0wl00 


11 


27 


2 


0- 


100 


mocr04 


36, 


726 


81 


803 


33 


446 


99.02 


gnlle41d2w0wl00 


11 


41 


2 


0- 


100 


mocr04 


120, 


393 


309 


961 


188 


619 


99.53 


gnl2ol6d2w0wl00 


12 


16 


2 


0- 


100 


mocr04 


1, 


750 




162 




161 


99.35 


gnl2o24d2w0wl00 


12 


24 


2 


0- 


100 


mocr04 


62, 


618 


92 


813 


47 


870 


96.49 


gnl2G33d2w0wl00 


12 


33 


2 


0- 


100 


mocr04 


78, 


539 


192 


122 


101 


544 


99.00 


gnl3el9d2w0wl00 


13 


19 


2 


0- 


100 


mocr04 


9, 


703 


3 


255 


2 


467 


94.59 


gnl3e29d2w0wl00 


13 


29 


2 


0- 


100 


mocr04 


68, 


776 


164 


617 


48 


053 


99.22 


gnl3o39d2w0wl00 


13 


39 


2 


0- 


100 


mocr04 


138, 


832 


315 


881 


197 


263 


99.18 
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Table 4.9. Calibration test set: Algorithm 88 (Boundary Calculation) with 2 criteria. All 
times under 1 second. 



Name 


Nodes 


Edges 




Weight 
range 


Computer 




^Spanning 
Trees 




:?^Projected 
Trees 


Boundary 
points 


gn9el8d2w0w20 


9 


18 





- 20 


moci'02 


2 


981 


1 


742 


20 


gn9e27d2w0w20 


9 


27 





- 20 


mocr02 


372 


320 


6 


346 


29 


gnl0e22d2w0w20 


10 


22 





- 20 


mocr02 


53 


357 


3 


957 


35 


gnl0e28d2w0w20 


10 


28 





- 20 


mocr02 


800 


948 


7 


131 


33 


gnl0e33d2w0w20 


10 


33 





- 20 


moci'02 


3,584 


016 


10 


833 


10 


gnllel3d2w0w20 


11 


13 





- 20 


mocr02 




41 




41 


10 


gnlle20d2w0w20 


11 


20 





- 20 


mocr02 


6 


939 


2 


173 


20 


gnlle27d2w0w20 


11 


27 





- 20 


mocr02 


284 


730 


7 


515 


27 


gnlle41d2w0w20 


11 


41 





- 20 


mocr02 


90, 922 


271 


16 


457 


38 


gnl2el6d2w0w20 


12 


16 





- 20 


mocr02 




162 




154 


13 


gnl2e24d2w0w20 


12 


24 





- 20 


mocr02 


208 


380 


5 


647 


25 


gnl2e33d2w0w20 


12 


33 





- 20 


mocr02 


4, 741 


624 


9 


364 


35 


gnl3el9d2w0w20 


13 


19 





- 20 


mocr02 


3 


401 


1 


608 


21 


gnl3e29d2w0w20 


13 


29 





- 20 


mocr02 


1, 543 


340 


8 


474 


36 


gnl3e39d2w0w20 


13 


39 





— 20 


mocr02 


131, 807 


934 


18 


468 


52 


gn9el8d2w0wl00 


9 


18 


- 


- 100 


mocr02 


2 


981 


2 


858 


19 


gn9e27d2w0wl00 


9 


27 


- 


~ 100 


mocr02 


372 


320 


89 


092 


31 


gnl0c22d2w0wl00 


10 


22 


- 


- 100 


mocr02 


53 


357 


37 


204 


26 


gnl0c28d2w0wl00 


10 


28 


- 


- 100 


mocr02 


800 


948 


101 


334 


31 


gnl0c33d2w0wl00 


10 


33 


0- 


- 100 


mocr02 


3,584 


016 


166 


427 


37 


gnllcl3d2w0wl00 


11 


13 


0- 


- 100 


mocr02 




41 




41 


9 


gnllc20d2w0wl00 


11 


20 


- 


- 100 


mocr02 


6 


939 


6 


580 


22 


gnllc27d2w0wl00 


11 


27 


- 


- 100 


mocr02 


284 


730 


81 


803 


30 


gnllc41d2w0wl00 


11 


41 


- 


~ 100 


mocr02 


90, 922 


271 


309 


961 


41 


gnl2cl6d2w0wl00 


12 


16 


- 


- 100 


mocr02 




162 




162 


11 


gnl2e24d2w0wl00 


12 


24 


- 


- 100 


mocr02 


208 


380 


92 


813 


29 


gnl2e33d2w0wl00 


12 


33 


- 


- 100 


mocr02 


4, 741 


624 


192 


122 


32 


gnl3el9d2w0wl00 


13 


19 


- 


- 100 


mocr02 


3 


401 


3 


255 


22 


gnl3c29d2w0wl00 


13 


29 


- 


- 100 


mocr02 


1, 543 


340 


164 


617 


37 


gnl3c39d2w0wl00 






- 


- lUU 


mocr02 


131, 807 




315 


ool 


A 


gn9ol8d2w0wl000 


9 


18 


- 


1000 


mocr02 


2 


981 


2 


981 


26 


gn9o27d2w0wl000 


9 


27 


- 


1000 


mocr02 


372 


320 


364 


382 


31 


gnl0e22d2w0wl000 


10 


22 


- 


1000 


mocr02 


53 


357 


52 


990 


37 


gnl0c28d2w0wl000 


10 


28 


- 


1000 


mocr02 


800 


948 


756 


013 


9 


gnl0o33d2w0wl000 


10 


33 


0- 


1000 


mocr02 


3,584 


016 


2,726 


287 


21 


gnllcl3d2w0wl000 


11 


13 


0- 


1000 


mocr02 




41 




41 


28 


gnllo20d2w0wl000 


11 


20 


0- 


1000 


mocr02 


6 


939 


6 


921 


42 


gnllc27d2w0wl000 


11 


27 


0- 


1000 


mocr02 


284 


730 


279 


308 


13 


gnllc41d2w0wl000 


11 


41 


0- 


1000 


mocr02 


90, 922 


271 


13,884 


793 


31 


gnl2cl6d2w0wl000 


12 


16 


0- 


1000 


moci'02 




162 




162 


32 


gnl2c24d2w0wl000 


12 


24 


0- 


1000 


mocr02 


208 


380 


205 


690 


18 


gnl2o33d2w0wl000 


12 


33 


0- 


1000 


mocr02 


4,741 


624 


3,680 


313 


36 


gnl3cl9d2w0wl000 


13 


19 


0- 


1000 


mocr02 


3 


401 


3 


401 


43 


gnl3c29d2w0wl000 


13 


29 


0- 


1000 


mocr02 


1,543 


340 


1,396 


180 


17 


gnl3e39d2w0wl000 


13 


39 


0- 


1000 


mocr02 


131,807 


934 


15,037 


589 


26 
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Table 4.10. Calibration test set; Heuristic 91 (Boundary and Triangular Region Pareto Test) 
with 2 criteria. Internally, Heuristic 86 (PT) used 10 searches per point, 100 pivot limit. 



Name 


Nodes 


Edges 




Weights 


Computer 


Seconds 


# Pareto Optimum 


#Computed 
Pareto Optimum 




Percent 


giiyc ±oc_izwuwiu 


9 


18 


- 


- 20 


uzzy 


40 


15 


15 


100 


00 


giiycz 1 Qzwuwzu 


Q 


97 


u - 




uzzy 


11 


1 9 


1 9 


lUU 


00 

yjyj 


(vn 1 rinOOi-nnrriiiii-on 

glliUCZZQZWUWZU 


1 n 


22 


U - 


- 90 


Fuzzy 




1 9 


1 9 


lUU 


00 

yjyj 


gn iUCZoQ^WUWZU 


1 n 




U - 


~ 90 


uzzy 


q 


1 


1 


lUU 


00 

yjyj 


gn iucooQzwuwzu 


1 n 

-L yj 




U - 


- 90 


uzzy 


59 


21 


21 


lUU 


00 

yjyj 


gni iciouzwuwzu 


1 1 
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Table 4.12. Calibration test set: Local Search Heuristic 82 minimizing (x — x)^, where for 
each instance, x is an interior integer point that is the image of some base. Here d = 2 and all 
tests run on Fuzzy. Shown are the number of spanning trees and projected bases of each instance 
and the number of successes out of 1000 Local Searches and seconds to perform the 1000 searches. 
We also repeat the exact number of spanning trees and exact number of projected spanning trees 
and the proportion (given as percentage) between them. 
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Table 4.13. Calibration test set: Local Search Heuristic 82 minimizing (x — x)^, where for 
each instance, x is an interior integer point that is the image of some base. Here d = 3 and all 
tests run on Fuzzy. Shown are the number of spanning trees and projected bases of each instance 
and the number of successes out of 1000 Local Searches and seconds to perform the 1000 searches. 
We also repeat the exact number of spanning trees and exact number of projected spanning trees 
and the proportion (given as percentage) between them. 
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Table A. 1.4:. Calibration test set: Local Search Heuristic 82 minimizing (x — x)^, where for 
each instance, x is an interior integer point that is the image of some base. Here d = 5 and all 
tests run on Fuzzy. Shown are the number of spanning trees and projected bases of each instance 
and the number of successes out of 1000 Local Searches and seconds to perform the 1000 searches. 
We also repeat the exact number of spanning trees and exact number of projected spanning trees 
and the proportion (given as percentage) between them. 
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Table 4.15. Calibration test set: Local Search Heuristic 82 minimizing (x — x)*, where for 
each instance, x is an interior integer point that is the image of some base. Here d = 2 and all 
tests run on Fuzzy. Shown are the number of spanning trees and projected bases of each instance 
and the number of successes out of 1000 Local Searches and seconds to perform the 1000 searches. 
We also repeat the exact number of spanning trees and exact number of projected spanning trees 
and the proportion (given as percentage) between them. 
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Table 4.16. Calibration test set: Local Search Heuristic 82 minimizing (x — x)*, where for 
each instance, x is an interior integer point that is the image of some base. Here d = 3 and all 
tests run on Fuzzy. Shown are the number of spanning trees and projected bases of each instance 
and the number of successes out of 1000 Local Searches and seconds to perform the 1000 searches. 
We also repeat the exact number of spanning trees and exact number of projected spanning trees 
and the proportion (given as percentage) between them. 
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Table 4.17. Calibration test set: Local Search Heuristic 82 minimizing (x — x)*, where for 
each instance, x is an interior integer point that is the image of some base. Here d = 5 and all 
tests run on Fuzzy. Shown are the number of spanning trees and projected bases of each instance 
and the number of successes out of 1000 Local Searches and seconds to perform the 1000 searches. 
We also repeat the exact number of spanning trees and exact number of projected spanning trees 
and the proportion (given as percentage) between them. 
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Table 4.24. Calibration test set: Local Search Heuristic 82 minimizing (x — x)^, where for 
each instance, x is an interior rational non-integer point. Here d = 2 and all tests run on Fuzzy. 
Shown are the number of spanning trees and projected bases of each instance and the number of 
successes out of 1000 Local Searches and seconds to perform the 1000 searches. We also repeat the 
exact number of spanning trees and exact number of projected spanning trees and the proportion 
(given as percentage) between them. 



Name 




^Spanning 
Trees 




#Projected 
Trees 


Percent 


^Successes 
out of 
1000 


Seconds 


gn9el8d2w0w20 


2 


981 


1 


742 


58.44 


144 


1 


gn9e27d2w0w20 


372 


320 


6 


346 


01.70 


273 


1 


gnl0e22d2w0w20 


53 


357 


3 


957 


07.42 


203 


2 


gnl0e28d2w0w20 


800 


948 


7 


131 


00.89 


364 


2 


gnl0e33d2w0w20 


3,584 


016 


10 


833 


00.30 


216 


2 


gnllel3d2w0w20 




41 




41 


100.00 


252 


1 


gnlle20d2w0w20 


6 


939 


2 


173 


31.32 


125 


2 


gnlle27d2w0w20 


284 


730 


7 


515 


02.64 


181 


1 


gnlle41d2w0w20 


90,922 


271 


16 


457 


00.02 


318 


2 


gnl2el6d2w0w20 




162 




154 


95.06 


240 


1 


gnl2e24d2w0w20 


208 


380 


5 


647 


02.71 


166 


2 


gnl2e33d2w0w20 


4,741 


624 


9 


364 


00.20 


219 


2 


gnl3el9d2w0w20 


3 


401 


1 


608 


47.28 


78 


1 


gnl3e29d2w0w20 


1,543 


340 


8 


474 


00.55 


199 


2 


gnl3e39d2w0w20 


131,807 


934 


18 


468 


00.01 


232 


3 


gn9el8d2w0wl00 


2 


981 


2 


858 


95.87 


440 


1 


gn9e27d2w0wl00 


372 


320 


89 


092 


23.93 





1 


gnl0e22d2w0wl00 


53 


357 


37 


204 


69.73 


47 


2 


gnl0e28d2w0wl00 


800 


948 


101 


334 


12.65 





2 


gnl0e33d2w0wl00 


3,584 


016 


166 


427 


04.64 


16 


2 


gnllel3d2w0wl00 




41 




41 


100.00 


283 


1 


gnlle20d2w0wl00 


6 


939 


6 


580 


94.83 





1 


gnlle27d2w0wl00 


284 


730 


81 


803 


28.73 





2 


gnlle41d2w0wl00 


90,922 


271 


309 


961 


00.34 


21 


3 


gnl2el6d2w0wl00 




162 




162 


100.00 


526 


1 


gnl2e24d2w0wl00 


208 


380 


92 


813 


44.54 





2 


gnl2e33d2w0wl00 


4,741 


624 


192 


122 


04.05 


7 


2 


gnl3el9d2w0wl00 


3 


401 


3 


255 


95.71 


78 


1 


gnl3e29d2w0wl00 


1,543 


340 


164 


617 


10.67 


8 


2 


gnl3e39d2w0wl00 


131,807 


934 


315 


881 


00.24 


34 


3 


gn9el8d2w0wl000 


2 


981 


2 


981 


100.00 


75 


1 


gn9e27d2w0wl000 


372 


320 


364 


382 


97.87 





1 


gnl0e22d2w0wl000 


53 


357 


52 


990 


99.31 





1 


gnl0e28d2w0wl000 


800 


948 


756 


013 


94.39 





2 


gnl0e33d2w0wl000 


3,584 


016 


2,726 


287 


76.07 





2 


gnllel3d2w0wl000 




41 




41 


100.00 


766 





gnlle20d2w0wl000 


6 


939 


6 


921 


99.74 


12 


1 


gnlle27d2w0wl000 


284 


730 


279 


308 


98.10 





1 


gnlle41d2w0wl000 


90,922 


271 


13,884 


793 


15.27 





2 


gnl2el6d2w0wl000 




162 




162 


100.00 


355 


1 


gnl2e24d2w0wl000 


208 


380 


205 


690 


98.71 


2 


2 


gnl2e33d2w0wl000 


4,741 


624 


3,680 


313 


77.62 





3 


gnl3el9d2w0wl000 


3 


401 


3 


401 


100.00 


1 


1 


gnl3e29d2w0wl000 


1,543 


340 


1,396 


180 


90.46 


1 


2 


gnl3e39d2w0wl000 


131,807 


934 


15,037 


589 


11.41 





3 
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Table 4.25. Calibration test set: Local Search Heuristic 82 minimizing (x — x)^, where for 
each instance, x is an interior rational non-integer point. Here d = 3 and all tests run on Fuzzy. 
Shown are the number of spanning trees and projected bases of each instance and the number of 
successes out of 1000 Local Searches and seconds to perform the 1000 searches. We also repeat the 
exact number of spanning trees and exact number of projected spanning trees and the proportion 
(given as percentage) between them. 



Name 








u 

O M 
0) 

OJ 

* 


Percent 


T^Succesf 
out of 
1000 


Seconds 


gn9el8d3w0w20 


2 


981 


2 


905 


97.45 


62 


1 


gn9e27d3w0w20 


372 


320 


134 


039 


36.00 





1 


gnl0e22d3w0w20 


53 


357 


42 


887 


80.38 


18 


2 


gnl0e28d3w0w20 


800 


948 


238 


529 


29.78 





2 


gnl0e33d3w0w20 


3, 584 


016 


411 


730 


11.49 





3 


gnllel3d3w0w20 




41 




41 


100.00 


183 





gnlle20d3w0w20 


6 


939 


6 


603 


95.16 


9 


1 


gnlle27d3w0w20 


284 


730 


146 


672 


51.51 


1 


2 


gnlle41d3w0w20 


90,922 


271 


959 


469 


01.06 


14 


2 


gnl2el6d3w0w20 




162 




162 


100.00 


206 


1 


gnl2e24d3w0w20 


208 


380 


111 


201 


53.36 


32 


2 


gnl2e33d3w0w20 


4, 741 


624 


470 


609 


09.93 


25 


2 


gnl3el9d3w0w20 


3 


401 


3 


358 


98.74 


17 


1 


gnl3e29d3w0w20 


1,543 


340 


264 


949 


17.18 


23 


2 


gnl3e39d3w0w20 


131,807 


934 


930 


322 


00.71 


14 


3 


gn9el8d3w0wl00 


2 


981 


2 


981 


100.00 


33 


1 


gn9e27d3w0wl00 


372 


320 


367 


313 


98.66 





2 


gnl0e22d3w0wl00 


53 


357 


53 


289 


99.87 


2 


1 


gnl0e28d3w0wl00 


800 


948 


786 


781 


98.23 





2 


gnl0e33d3w0wl00 


3,584 


016 


3,351 


096 


93.01 





2 


gnllel3d3w0wl00 




41 




41 


100.00 


700 





gnlle20d3w0wl00 


6 


939 


6 


939 


100.00 





2 


gnlle27d3w0wl00 


284 


730 


281 


303 


98.80 





1 


gnlle41d3w0wl00 


90,922 


271 


35, 943 


327 


39.53 





3 


gnl2el6d3w0wl00 




162 




162 


100.00 


825 


1 


gnl2e24d3w0wl00 


208 


380 


207 


143 


99.41 


3 


2 


gnizeooaowOwiUO 


4,741 


624 


4,740 


880 


99.98 





3 


gnl3el9d3w0wl00 


3 


401 


3 


401 


100.00 


59 


1 


gnl3e29d3w0wl00 


1,543 


340 


1,484 


719 


96.20 





3 


gnl3e39d3w0wl00 


131,807 


934 


44, 757 


592 


33.96 





4 


gn9el8d3w0wl000 


2 


981 


2 


981 


100.00 


235 


1 


gn9e27d3w0wl000 


372 


320 


372 


320 


100.00 





1 


gnl0e22d3w0wl000 


53 


357 


53 


357 


100.00 


10 


1 


gnl0e28d3w0wl000 


800 


948 


800 


946 


99.99 





2 


gnl0e33d3w0wl000 


3,584 


016 


3,583 


757 


99.99 





2 


gnllel3d3w0wl000 




41 




41 


100.00 


410 


1 


gnlle20d3w0wl000 


6 


939 


6 


939 


100.00 





1 


gnlle27d3w0wl000 


284 


730 


284 


730 


100.00 


1 


2 


gnlle41d3w0wl000 


90,922 


271 


90, 699 


181 


99.75 





3 


gnl2el6d3w0wl000 




162 




162 


100.00 


583 


1 


gnl2e24d3w0wl000 


208 


380 


208 


356 


99.99 





2 


gnl2e33d3w0wl000 


4,741 


624 


4,740 


880 


99.98 





3 


gnl3el9d3w0wl000 


3 


401 


3 


401 


100.00 


51 


1 


gnl3e29d3w0wl000 


1,543 


340 


1,543 


304 


99.99 





3 


gnl3e39d3w0wl000 


131,807 


934 


131,464 


478 


99.74 





3 
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Table 4.26. Calibration test set: Local Search Heuristic 82 minimizing (x — x)^, where for 
each instance, x is an interior rational non-integer point. Here d = 5 and all tests run on Fuzzy. 
Shown are the number of spanning trees and projected bases of each instance and the number of 
successes out of 1000 Local Searches and seconds to perform the 1000 searches. We also repeat the 
exact number of spanning trees and exact number of projected spanning trees and the proportion 
(given as percentage) between them. 



Name 




'3 

ft ^ 

*S 




u 

-r-5 Q 

o 
u u 

* 


Percent 


out of 
1000 


Seconds 


gn9el8d5w0wl 


2 


981 


1 


099 


36.87 


624 


2 


gn9e27d5w0wl 


372 


320 


8 


205 


2.20 


766 


2 


gnl0e22d5w0wl 


53 


357 


4 


746 


8.89 


693 


2 


gnl0e28d5w0wl 


800 


948 


10 


898 


1.36 


774 


2 


gnl0e33d5w0wl 


3, 584 


016 


15 


482 


0.43 


956 


3 


gnllelSdSwOwl 




41 




27 


65.85 


504 





gnlle20d5w0wl 


6 


939 


1 


866 


26.89 


529 


2 


gnlle27d5w0wl 


284 


730 


5 


828 


2.05 


608 


2 


gnlle41d5w0wl 


90,922 


271 


36 


847 


0.04 


953 


3 


gnl2el6d5w0wl 




162 




116 


71.61 


286 


1 


gnl2e24d5w0wl 


208 


380 


6 


738 


3.23 


860 


3 


gnl2e33d5w0wl 


4, 741 


624 


20 


857 


0.44 


957 


3 


gn 13e 19d5 wOwl 


3 


401 


1 


149 


33.78 


448 


1 


gnl3e29d5w0wl 


1, 543 


340 


111 


116 


7.20 


746 


3 


gn 1 3e39d5 wOwl 


131,807 


934 


34 


392 


0.03 


623 


5 


gn9el8d5w0w2 


2 


981 


2 


276 


76.35 


308 


1 


gn9e27d5w0w2 


372 


320 


43 


029 


11.55 


323 


2 


gnl0e22d5w0w2 


53 


357 


14 


623 


27.40 


97 


2 


gnl0e28d5w0w2 


800 


948 


66 


190 


8.26 


268 


2 


gnl0c33d5w0w2 


3, 584 


016 


105 


309 


2.93 


189 


3 


gnllol3d5w0w2 




41 




41 


100.00 


706 


1 


gnllo20d5w0w2 


6 


939 


3 


761 


54.20 


177 


1 


gnlle27d5w0w2 


284 


730 


39 


292 


13.79 


257 


3 


gnlle41d5w0w2 


90,922 


271 


290 


555 


0.31 


430 


4 


gnl2el6d5w0w2 




162 




160 


98.76 


245 


1 


gnl2e24d5w0w2 


208 


380 


34 


057 


16.34 


146 


2 


gn i 2eooao wUwz 


4,741 


624 


120 


211 


2.53 


404 


3 


gnl3el9d5w0w2 


3 


401 


2 


613 


76.83 


28 


2 


gnl3e29d5w0w2 


1, 543 


340 


98 


285 


6.36 


288 


4 


gnl3e39d5w0w2 


131,807 


934 


348 


703 


0.26 


361 


4 


gn9e 18d5 w0w5 


2 


981 


2 


960 


99.29 


32 


1 


gn9e27d5w0w5 


372 


320 


271 


048 


72.79 





2 


gnl0e22d5w0w5 


53 


357 


49 


463 


92.70 





1 


gnl0e28d5w0w5 


800 


948 


493 


565 


61.62 


9 


2 


gnl0e33d5w0w5 


3,584 


016 


1,294 


875 


36.12 


1 


3 


gnllcl3d5w0w5 




41 




41 


100.00 


574 





gnlle20d5w0w5 


6 


939 


5 


975 


86.10 


14 


1 


gnlle27d5w0w5 


284 


730 


222 


974 


78.31 





2 


gnlle41d5w0w5 


90,922 


271 


4,711 


354 


5.18 


33 


4 


gnl2el6d5w0w5 




162 




162 


100.00 


194 





gnl2e24d5w0w5 


208 


380 


177 


845 


85.34 





2 


gn 1 2e33d5 w0w5 


4,741 


624 


1,489 


751 


31.41 


27 


4 


gnl3el9d5w0w5 


3 


401 


3 


381 


99.41 


48 


1 


gnl3e29d5w0w5 


1,543 


340 


787 


196 


51.00 


25 


4 


gnl3e39d5w0w5 


131,807 


934 


7,737 


684 


5.87 


16 


5 
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Table 4.27. Calibration test set: Local Search Heuristic 82 minimizing (x — x)*, where for 
each instance, x is an interior rational non-integer point. Here d = 2 and all tests run on Fuzzy. 
Shown are the number of spanning trees and projected bases of each instance and the number of 
successes out of 1000 Local Searches and seconds to perform the 1000 searches. We also repeat the 
exact number of spanning trees and exact number of projected spanning trees and the proportion 
(given as percentage) between them. 



Name 




^Spanning 
Trees 




#Projected 
Trees 


Percent 


^Successes 
out of 
1000 


Seconds 


gn9el8d2w0w20 


2 


981 


1 


742 


58.44 


161 


1 


gn9e27d2w0w20 


372 


320 


6 


346 


01.70 


211 


1 


gnl0e22d2w0w20 


53 


357 


3 


957 


07.42 


194 


1 


gnl0e28d2w0w20 


800 


948 


7 


131 


00.89 


363 


1 


gnl0e33d2w0w20 


3,584 


016 


10 


833 


00.30 


254 


2 


gnllel3d2w0w20 




41 




41 


100.00 


210 


1 


gnlle20d2w0w20 


6 


939 


2 


173 


31.32 


124 


1 


gnlle27d2w0w20 


284 


730 


7 


515 


02.64 


234 


1 


gnlle41d2w0w20 


90,922 


271 


16 


457 


00.02 


390 


2 


gnl2el6d2w0w20 




162 




154 


95.06 


260 


1 


gnl2e24d2w0w20 


208 


380 


5 


647 


02.71 


166 


1 


gnl2e33d2w0w20 


4,741 


624 


9 


364 


00.20 


224 


2 


gnl3el9d2w0w20 


3 


401 


1 


608 


47.28 


82 


2 


gnl3e29d2w0w20 


1,543 


340 


8 


474 


00.55 


172 


1 


gnl3e39d2w0w20 


131,807 


934 


18 


468 


00.01 


240 


3 


gn9el8d2w0wl00 


2 


981 


2 


858 


95.87 


373 


1 


gn9e27d2w0wl00 


372 


320 


89 


092 


23.93 


5 


1 


gnl0e22d2w0wl00 


53 


357 


37 


204 


69.73 


53 


1 


gnl0e28d2w0wl00 


800 


948 


101 


334 


12.65 


27 


2 


gnl0e33d2w0wl00 


3,584 


016 


166 


427 


04.64 


12 


1 


gnllel3d2w0wl00 




41 




41 


100.00 


330 





gnlle20d2w0wl00 


6 


939 


6 


580 


94.83 





1 


gnlle27d2w0wl00 


284 


730 


81 


803 


28.73 





1 


gnlle41d2w0wl00 


90,922 


271 


309 


961 


00.34 


18 


3 


gnl2el6d2w0wl00 




162 




162 


100.00 


537 


1 


gnl2e24d2w0wl00 


208 


380 


92 


813 


44.54 





1 


gnl2e33d2w0wl00 


4,741 


624 


192 


122 


04.05 


5 


2 


gnl3el9d2w0wl00 


3 


401 


3 


255 


95.71 


89 


1 


gnl3e29d2w0wl00 


1,543 


340 


164 


617 


10.67 


18 


1 


gnl3e39d2w0wl00 


131,807 


934 


315 


881 


00.24 


34 


3 


gn9el8d2w0wl000 


2 


981 


2 


981 


100.00 


71 


1 


gn9e27d2w0wl000 


372 


320 


364 


382 


97.87 





2 


gnl0e22d2w0wl000 


53 


357 


52 


990 


99.31 





1 


gnl0e28d2w0wl000 


800 


948 


756 


013 


94.39 





2 


gnl0e33d2w0wl000 


3,584 


016 


2,726 


287 


76.07 





1 


gnllel3d2w0wl000 




41 




41 


100.00 


434 





gnlle20d2w0wl000 


6 


939 


6 


921 


99.74 


12 


1 


gnlle27d2w0wl000 


284 


730 


279 


308 


98.10 





2 


gnlle41d2w0wl000 


90,922 


271 


13,884 


793 


15.27 





2 


gnl2el6d2w0wl000 




162 




162 


100.00 


346 


1 


gnl2e24d2w0wl000 


208 


380 


205 


690 


98.71 





1 


gnl2e33d2w0wl000 


4,741 


624 


3,680 


313 


77.62 





2 


gnl3el9d2w0wl000 


3 


401 


3 


401 


100.00 


5 


1 


gnl3e29d2w0wl000 


1,543 


340 


1,396 


180 


90.46 





1 


gnl3e39d2w0wl000 


131,807 


934 


15,037 


589 


11.41 





3 
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Table 4.28. Calibration test set: Local Search Heuristic 82 minimizing (x — x)*, where for 
each instance, x is an interior rational non-integer point. Here d = 3 and all tests run on Fuzzy. 
Shown are the number of spanning trees and projected bases of each instance and the number of 
successes out of 1000 Local Searches and seconds to perform the 1000 searches. We also repeat the 
exact number of spanning trees and exact number of projected spanning trees and the proportion 
(given as percentage) between them. 



a 



Name 








u 

o m 
0) 

OJ 
u u 

* 


Percent 


T^Succesf 
out of 
1000 


Seconds 


gn9el8d3w0w20 


2 


981 


2 


905 


97.45 





1 


gn9e27d3w0w20 


372 


320 


134 


039 


36.00 


8 


2 


gnl0e22d3w0w20 
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Table 4.29. Calibration test set: Local Search Heuristic 82 minimizing (x — x)*, where for 
each instance, x is an interior rational non-integer point. Here d = 5 and all tests run on Fuzzy. 
Shown are the number of spanning trees and projected bases of each instance and the number of 
successes out of 1000 Local Searches and seconds to perform the 1000 searches. We also repeat the 
exact number of spanning trees and exact number of projected spanning trees and the proportion 
(given as percentage) between them. 
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Table 4.36. Two experimental design problems P20 and P28 and comparing 
run times of IBM Blue Gene/L vs Fuzzy to find all projected bases. We used 
Heuristic 92 (DFBFS) with 100 searches, boundary retry limit 100, interior 
retry limit 10, 000 and truncation depth 8. 



Name 




P2Q 




-P28 


Loiumns(7y:matroid eiementsj 




100 




100 


Rows(matroid rank) 




20 




28 


Projected bases 


13 


816 


19, 


193 


Seconds for IBM Blue Gene/L 


23 


773 


56, 


763 


^ Projected bases found by run 1 of Heuristic 92 


13 


757 


19, 


193 


Seconds for run 1 of Heuristic 92 




923 


2, 


557 


^ Projected bases found by run 2 of Heuristic 92 


13 


816 


19, 


191 


Seconds for run 2 of Heuristic 92 




945 


2, 


645 



4.3. COMPUTATIONAL RESULTS 



124 



Table 4.37. Experimental design problem P20 ^-nd comparing run times 
of IBM Blue Gene/L vs Fuzzy to find all projected bases. We partitioned 
the region of feasible projected bases into 12 disjoint regions and ran 12 
separate concurrently running instances of Heuristic 86 (PT) on Fuzzy to 
find as many projected bases as possible. When running Heuristic 86 (PT) 
we used Heuristic 85 (TS) and tried each point 10 times with a Tabu limit 
of 20. 



Region 




)jected 
5 Found 






Basel 






ID 
Ul 


Boxl 


1, 


044 


9,544 


Box2 


1, 


122 


5,086 


Box3 


1, 


122 


4,966 


Box4 


1, 


135 


8,066 


Box5 


1, 


122 


5,414 


Box6 


1, 


122 


4,540 


Box7 


1, 


122 


4,528 


Box8 


1, 


188 


5,614 


Box9 


1, 


147 


9,181 


BoxlO 


1, 


124 


5,232 


Boxll 


1, 


124 


5,252 


Boxl2 


1, 


249 


8,173 


TOTAL 


13, 


816 


(Longest time)9,544 
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Table 4.38. Machines Used in Experiments. 





#CPUS 


Clock Soeed 


A/fpTnorv 


Fuzzy- 


16 AMD 


2400 MHz 


64GB 


Truth 


16 AMD 


2400 MHz 


64GB 


MathocrOl 


2 AMD 


2000 MHz 


4GB 


Mathocr02 


2 AMD 


2000 MHz 


4GB 


Mathocr03 


2 AMD 


2000 MHz 


4GB 


Mathocr04 


2 AMD 


2000 MHz 


12GB 
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